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ABSTRACT 

Deep long-slit optical spectrophotometric observations are presented for 25 Galactic bulge 
planetary nebulae (GBPNe) and 6 Galactic disk planetary nebulae (GDPNe). The spectra, 
i c| combined with archival ultraviolet spectra obtained with the International Ultraviolet Ex- 

Q_C plorer (IUE) and infrared spectra obtained with the Infrared Space Observatory (ISO), have 

JL | been used to carry out a detailed plasma diagnostic and element abundance analysis utilizing 

both collisional excited lines (CELs) and optical recombination lines (ORLs). 

■ Comparisons of plasma diagnostic and abundance analysis results obtained from CELs 
| and from ORLs reproduce many of the patterns previously found for GDPNe. In particular we 

'— 1 . show that the large discrepancies between electron temperatures (T e 's) derived from CELs and 

from ORLs appear to be mainly caused by abnormally low values yielded by recombination 
lines and/or continua. Similarly, the large discrepancies between heavy element abundances 
deduced from ORLs and from CELs are largely caused by abnormally high values obtained 
from ORLs, up to tens of solar in extreme cases. It appears that whatever mechanisms are 
ly-j , causing the ubiquitous dichotomy between CELs and ORLs, their main effects are to enhance 

■ the emission of ORLs, but hardly affect that of CELs. It seems that heavy element abundances 

deduced from ORLs may not reflect the bulk composition of the nebula. Rather, our analysis 
suggests that ORLs of heavy element ions mainly originate from a previously unseen compo- 
nent of plasma of T e 's of just a few hundred Kelvin, which is too cool to excite any optical 
and UV CELs. 

• • . We find that GBPNe are on the average 0. 1-0.2 dex more metal-rich than GDPNe but 

. ' have a mean C/O ratio that is approximately 0.2 dex lower. By comparing the observed rela- 

k> , tive abundances of heavy elements with recent theoretical predictions, we show that GBPNe 

5_j ■ probably evolved from a relatively metal-rich environment of initial Z ~ 0.013, compared to 

an initial Z < 0.008 for GDPNe. In addition, we find that GBPNe tend to have more massive 
progenitor stars than GDPNe. GBPNe are found to have an average magnesium abundance 
about 0.13 dex higher than GDPNe. The latter have a mean magnesium abundance almost 
identical to the solar value. The enhancement of magnesium in GBPNe and the large [a/Fe] 
ratios of bulge giants suggest that the primary enrichment process in the bulge was Type II 
SNe. PN observations yield a Ne/O abundance ratio much higher than the solar value, sug- 
gesting that the solar neon abundance may have been underestimated by 0.2 dex. 
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1 INTRODUCTION 

Knowledge of the chemical composition of the Galactic bulge 
is of paramount importance for the understanding of the history 
of fo r mation and evo lution of the Galaxy (e.g. Ivan den Berghl 
ll99dlMcWilliamlll997l) . Wh ile early work sug gested that bulge 
stars are super metal- rich (eg. Iwhitford & Richll 19831 : iRichll 19881 : 
lGeisler&Frielll992h . more recent high resolution spectroscopic 



studies showed that the bulge i s actually slightly iron-po or com- 
pared to the solar neig hborhood I McWillia m~& Richl 19941) . a result 
supported by the analysis of low resolution spectroscopic and pho- 
to metric observations for a large sample of bulge K and M giants 
bv llbata & Gilmorej < 1995b . In addition, a careful metallicity anal- 
ysis of low resolution spectr a of integrat ed light from the Galactic 
bulge in Baade's window bv lldiart et alj i 19961) suggested that the 
bulge is enhanced in a-elements (represented by magnesium). 
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All the aforementioned studies were based on analysis of light 
from bulge stars and yielded essentially no information of abun- 
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dances of light elements such as carbon and oxygen. In this con- 
text, abundance analyses using GBPNe as probes can play an im- 
portant role. Planetary nebulae (PNe) are products of low- and 
intermediate-mass stars (LIMS). Ionizing stars of PNe have typi- 
cal luminosities of several hundreds to thousands solar luminosity. 
Much of the energy is emitted in a few strong narrow emission 
lines, making them easily observable to large distances. In the op- 
tical and infrared, many PNe in the Galactic bulge can be readily 
observed even though the bulge is heavily obscured by interven- 
ing dust grains in the Galactic disk. For most elements, PNe are 
believed to preserve much of the original composition of the in- 
terstellar medium (ISM) at the time of their formation. Elemental 
abundances derived from PN observations can therefore be com- 
pared and related to results from stellar abundance analyses. Pos- 
sible exceptions are light elements such as He, N and C, whose 
abundances can potentially be modified by various dredge-up pro- 
cesses, depending on mass and initial metallicity of the progenitor 
star. Measurements of abundances of those elements in PNe thus 
shed light on the nucleosynthesis and dredge processes that oc- 
curred during the late evolutionary stages of their progenitor stars. 

Hitherto, spectroscopic abundance determinations have been 
carried out for approximately 250 GBPNe. The analyses were al- 
most invariably based on observations of rela t ively strong CEL s 
dRatag et al.|[l992l 1 19971 ; lEscudero et al]|2004 lExter et alHioj) . 
The results show that, albeit with large scatter, GBPNe have av- 
erage abundances and abundance distributions similar to GDPNe 
for elements unaffected by evolution of their progenitor stars, such 
as oxygen and other heavier elements. For helium and nitrogen, 
whose abundances may have been altered by the various dredge-up 
processes, studies so far have yielded different results. Part of the 
discrepancies may be due to differences in PN populations sampled 
by individual studies. 

The heavy extinction towards the Galactic bulge prohibits ex- 
tensive ultraviolet (UV) spectroscopy of GBPNe. Only a very re- 
strictive number of UV spectra, obtained mainly with the IUE, are 
available for GBPNe, mostly located within Baade's window. This 
poses a serious limitation for the study of some key elements, such 
as carbon for which the only observable strong CELs, the C IV 
AA1548J550 and C m] AA1907,1909 lines, fall in the ultraviolet 
(the infrared [C II] 158/^m line arises mainly from the photodisso- 
ciation region surrounding the ionized region, rather than from the 
latter itself). As a consequence, very little is known at the moment 
about the abundance and its distribution of this key element for the 
bulge population. Even with suitable UV observations, the analysis 
is complicated by the uncertain extinction law towards the Galactic 
bulge, which is known to deviate from that of the general ISM of 
the Galactic disk (Walton, Barlow & Clegg 1993; Liu et al. 2001). 

Advances in the past decade in observational techniques and 
in calculations of the basic atomic data, such as effective recombi- 
nation coefficients for non-hydrogenic ions with multiple valence 
electrons, have made it possible to determine nebular abundances 
using faint optical recombination lines from heavy element ions 
(e.g. Liu et al. 1995, 2000, 2001). Ionic abundances deduced from 
ORLs have the advantage that they are almost independent of the 
physical conditions, such as T e and electron density (N c ), of the 
nebula under study. Dominant ionic species of carbon typically 
found in photoionized gaseous nebulae possess a number of rela- 
tively bright ORLs, such as C II A4267, C III A4187 and A4650 and 
C IV 4658 (blended with the [Fe III] A4658 line in some objects), 
allowing measurements of ionic abundances of C 2+ , C 3+ and C 4+ . 
These lines are routinely detected in spectra of reasonable quality 
for bright Galactic PNe (e.g. Liu 1998). 



Deep, medium resolution optical spectroscopy, allowing de- 
terminations from ORLs of abundances of carbon, nitrogen, oxygen 
and neon, have been carried out for several dozen GDPNe (Tsamis 
et al. 2003, 2004; Liu et al. 2004a,b; Wesson, Liu & Barlow 2005). 
It is found that ionic abundances derived from ORLs are ubiqui- 
tously higher than values derived from CELs, i.e. abundance dis- 
crepancy factor (adf) > 1 for all the nebulae analyzed so far. The 
adf is found to vary from nebula to nebula and typically falls in 
the range 1.6-3.0, but with a significant tail extending to much 
higher values, reaching 70 in the most extreme case found so far 
(c.f. Liu 2003, 2006b for recent reviews). On the other hand, analy- 
ses also show that for a given nebula, adf 's derived for all abundant 
second-row elements, C, N, O and Ne, are of similar magnitude; in 
other words, heavy element abundance ratios, such as C/O, N/O and 
Ne/O are unaffected, regardless of the actual value of the adf (Liu 
et al. 1995, 2000, 2001; Mathis et al. 1998; Luo et al. 2001; Tsamis 
et al. 2004; Liu et al. 2004b). In stark contrast to those abundant 
second-row elements, ORL abundance determinations for magne- 
sium, the only third-row element that has been studied using an 
ORL, yield nearly a c onstant Mg/H abundance ratio almost identi- 
cal to the solar value tearlow et al J 20031 ; iLitill 2006bh . 

Detailed plasma diagnostics using ORLs indicate that ORLs 
may originate in regions of T c and iV c quite different from those 
where CELs are emi tted. While new obse rvations have confirmed 
the earlier results of lPeimbertl dl97ll) and lLiu & Danzigej dl993h 
that T c 's, derived from the Balmer discontinuity of the hydrogen 
recombination spectrum are systematically lower than values de- 
duced from the collisionally excited [O III] nebular to auroral for- 
bidden line ratio, careful analyses of the relative intensities of He I 
ORLs suggest that He I ORLs may arise from plasmas of T c 's even 
lower than indicated by the H I Balmer discontinuity. Much more 
importantly, careful comparisons of high quality measurements of 
observed relative intensities of O II ORLs in spectra of PNe exhibit- 
ing particularly large adf 's (> 5) clearly show that the lines arise 
from plasmas of T c ~ 1000 K, about an order of magnitude lower 
than values deduced from the [O III] forbidden line ratios. The ob- 
servations thus point to the presence of another component of cold 
plasma, previously unknown, embedded in the diffuse nebula. The 
gas has much higher metallicity and, because of the much enhanced 
cooling, a much lower T c . In addition, because of the very low T c , 
much lower than the typical excitation energies of optical and UV 
CELs, it emits essentially no optical and UV CELs, and is thus in- 
visible via those lines. The existence of such a cold, H-deficient 
component of plasma thus provides a natural explanation for the 
systematic higher heavy elemental abundances and lower T c 's de- 
duced from ORLs compared to those derived from CELs (Liu et 
al. 2000, 2006). The nature and origin of the high-metallicity gas 
remain elusive. Accurate determinations of the prevailing physical 
conditions (T c and iV c ) and elemental abundances will be essential. 

In this paper, we present deep, high-quality optical spectra for 
a sample of 25 GBPNe plus 6 GDPNe. The data are combined with 
archival IR and UV spectra to study the nebular physical conditions 
and elemental abundances using both CELs and ORLs. The sample, 
together with those published previously, mainly for GDPNe, bring 
the total number of PNe for which ORL elemental abundances have 
been determined to nearly a hundred. The purposes of the study 
are two folds: 1) to determine and characterize the distribution of 
adf's for a significant sample of GBPNe and compare it to that 
for GDPNe; 2) to determine carbon abundances and C/O ratios for 
the first time for a significant sample of GBPNe which are poorly 
known for the bulge population. The paper is organized as follows. 
In Section 2, we describe briefly our new optical observations and 



archival IR and UV data. Plasma diagnostics, ionic and elemental 
abundance determinations using CELs and ORLs are presented re- 
spectively in Sections 3 and 4. In Section 5, we compare our results 
with those in the literature. In Section 6, we present a statistical 
analysis of adf 's and other nebular physical properties. In Section 
7, we contrast the abundance patterns deduced for GDPNe and for 
GBPNe and compare them to theoretical models. We conclude by 
summarizing the main results in the last Section. 



2 TARGETS AND OBSERVATIONS 

In total 31 PNe in the direction of the Galactic center were ob- 
served. The basic parameters of the targets, including observed 
absolute H/3 flux, radio flux density at 5.0 GHz, angular diameter 
(value measured in the optical is adopted if available, otherwise 
radio measurement is used), heliocentric radial velocity and mag- 
nitude of the central star are given in Table Q] We assume that a 
PN probably belongs to the bulge population if it falls within 20 
degrees of the Galactic center, has an angular diameter less than 
20", a 5 GHz radio flux density of less than 50 mjy and a rela- 
tively large heliocentric radial velocity (Ratag etalJI 19971) . 25 of 
our targets satisfy these criteria. For the remaining 6 nebulae, H 1- 
35, M 1-29, M 1-61 and NGC 6567 have a radio flux density in 
excess of 70 mjy, and IC 4846 and He 2-1 18 fall nearly 30 degrees 
or more from the Galactic center. They are thus more likely belong- 
ing to the disk rather than the bulge population. We include these 
six PNe in our sample of GDPNe (c.f. Section[7}. 



2.1 Optical observations 

The spectra were obtained in 1995, 1996 and 2001 with the Eu- 
ropean Southern Observatory (ESO) 1.52 m telescope using the 
Boiler & Chivens (B&C) long-slit spectrograph. A journal of ob- 
servations is presented in Table [2] In 1995, the spectrograph was 
equipped with a Ford 2048x2048 15 /im x 15 CCD, which 
was superseded in 1996 by a UV-enhanced Loral 2048 x 2048 
15 /im x 15 ,um chip and in 2001 by a Loral 2688 x 2688 
15 /im x 15 fim chip. 

Two grating setups were used. The low resolution setting 
yielded a FWHM of ~ 4.5 or 6.0 A and covered the wavelength 
range from 3520 A to 7420 A. The high resolution setting gave a 
full width half maximum (FWHM) of ~ 0.9 or 1.5 A and covered 
the wavelength range 3994 - 4983 or 3500 - 4800 A. Typical ex- 
posure times were about 5 min for the low resolution setting. For 
the high spectral resolution, integration times of 20 to 30 min were 
used, and in most cases two or more frames were taken to increase 
the signal-to-noise (S/N) ratio and to facilitate removal of cosmic 
rays. The observations were aimed to detect at least the strongest 
ORLs of O II, i.e. the A4649 line from Multiplet V 1. In order to 
obtain strengths of the very strong emission lines, such as [O III] 
AA4959, 5007 and Ha which saturated on long exposures, short 
exposures of duration of 10 to 60 seconds were also obtained. A 
slit width of 2" was used for all nebular observations. However, in 
order to obtain the total flux of H/3, a short exposure was taken us- 
ing an 8" wide slit for each nebula for the low spectral resolution 
setup. 

All nebulae were observed with a fixed position long-slit, nor- 
mally passing through the nebular center and sampling the bright- 
est parts of the nebula. The only exception was Cn 1-5 for which 
scanned spectra were obtained by uniformly driving a long slit 
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Table 2. Journal of ESO 1.52-m observations. 



Name 


Date 
(UT) 


A-range 

(A) 


FWHM 

(A) 


Exp. Time 

(sec) 


Cn 1-5 


07/1996 


3520-7420 


4.5 


30,300 




07/1995 


3994^4983 


1.5 


2x1200 


Cn 2-1 


07/1995 


3520-7420 


4.5 


15,300 




07/1995 


3994^1983 


1.0 


2x1800 




06/2001 


3500-4805 


1.5 


1200 


H 1-35 


07/1995 


3520-7420 


4.5 


30,300 




07/1996 


3994^1983 


1.0 


2x1800 




06/2001 


3500-4805 


1.5 


1200 




06/2001 


4660-7260 


3.0 


900 


H 1-41 


07/1995 


3520-7420 


4.5 


30,2x300 




07/1996 


3994-4983 


1.5 


2x1200 


H 1-42 


07/1995 


3520-7420 


4.5 


30,300 




07/1996 


3994-4983 


1.0 


2x1800 


H 1-50 


07/1995 


3520-7420 


4.5 


30,300 




07/1995 


3994-4983 


1.0 


2x1800 




06/2001 


3500-4805 


1.5 


1200 


H 1-54 


07/1996 


3527-7431 


6.0 


60,300 




07/1996 


3994-4983 


1.5 


2x1800 




06/2001 


4660-7260 


1.5 


900 


He 2-118 


07/1996 


3520-7420 


4.5 


30,300 




07/1995 


3994^1983 


1.5 


2x1200 




06/2001 


3500-4805 


1.5 


1200 


IC 4699 


07/1995 


3520-7420 


4.5 


90,300 




07/1995 


3993^1979 


0.9 


2x1800 


IC 4846 


07/1995 


3520-7420 


4.5 


30,300 




07/1995 


3993^1980 


0.9 


2x1800 


M 1-20 


07/1996 


3527-7431 


6.0 


60,300 




07/1996 


3994^984 


1.5 


2x1800 




06/2001 


3500-4805 


1.5 


1800 


M 1-29 


07/1996 


3527-7431 


6.0 


60,600 


M 1-61 


07/1995 


3520-7420 


4.5 


20,300 




07/1995 


3993^1980 


0.9 


300,2x1800 


M2-4 


07/1995 


3520-7420 


4.5 


30,300 




07/1996 


3994^4984 


1.5 


2x1800 




06/2001 


3500-4805 


1.5 


1800 


M2-6 


07/1996 


3520-7420 


6.0 


30,300 




07/1996 


3994-^984 


1.5 


1800 


M2-23 


07/1995 


3520-7420 


4.5 


15,300 




07/1996 


3994-^984 


1.5 


2x1500 


M2-27 


07/1996 


3527-7431 


6.0 


60,600 


M2-31 


07/1996 


3527-7431 


6.0 


60,300 


M2-33 


07/1996 


3527-7431 


6.0 


60,300 




07/1996 


3994-^984 


1.5 


2x1800 


M2-39 


07/1995 


3520-7420 


4.5 


120,300 




07/1996 


3994-^984 


1.5 


2x1800 


M2-42 


07/1996 


3520-7431 


6.0 


60,300 


M3-7 


07/1996 


3527-7431 


6.0 


60,300 


M3-21 


07/1995 


3520-7420 


4.5 


10,300 




07/1996 


3994-4984 


0.9 


2x1800 




Uo/ZUUl 


jjUU— 


1 .J 


zxl IXjV) 


M3-29 


07/1996 


3520-7431 


6.0 


60,300 


M 3-32 


07/1995 


3530-7430 


4.5 


60,2x300 




07/1996 


3030-4050 


1.5 


2x1800 




07/1996 


3990-5000 


1.5 


2x1800 


M 3-33 


07/1995 


3520-7420 


4.5 


60,300 




07/1996 


3994-4984 


1.5 


2x1800 


NGC 6439 


07/1995 


3520-7420 


4.5 


30,300 




07/1996 


3994-4983 


1.5 


2x1200 




06/2001 


3500-4805 


1.5 


1200,2x1800 


NGC 6565 


07/1995 


3520-7420 


4.5 


90,300 




07/1995 


3993-4979 


0.9 


2x1800 
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Table 1. Basic data. Angular diameters are optical measurements except as specified. 



Name 


PNG 


Diam 


-logF(H/3) 


F(5GHz) 


RV 


B 


V 






(") 


(erg cm" 


~ 2 s" 1 ) 


(mjy ) 


(kms" 1 ) 


(mag) 


(mag) 






A92 


CKS92 


WL07 




D98 


T91.T89 




















Cn 1-5 


on? 9—00 4 


7 


11.25 


11.40 


44 


—29 


15.50 


15.20 


Cn 2-1 


9—04. 4 


2 4 


11.64 


11.63 


49 


—271 


* 




H 1-41 


J JO. / UH-.o 


Q 

7.D 


11.7° 


11.90 


12 


76 


16.30 


16.20 


H 1-42 


K7 9 04 S 


5 8 


11.7 a 


11.68 


40 


7Q 


17.30 




H 1-50 


JJO. / UJ.Z 


1 4 & 


11.68 


11.70 


J 1 .U 




* 


17.10 


H 1-54 


009 1 04 9 


4 8 


11.89 


11.87 


Jl.U 


1 1 u 


15.70 


15.40 


IC 4699 


148 0-11 8 


5.0 


11.69 


11.70 


20 


— 123 


14.82 


15.10 


M 1-20 


00 A 1 4-0 8 1 
UUU. 1 tuo . J 


1.7 


11.93 


11.94 




Q 1 
71 


17.70 


17.10 


M 2-4 


"?4Q 84-04 4 


5 


11.84 


11.94 


32 


— 184 


17.60 


17.00 


M 2-6 


JJJ.JtUO.J 


8 


12.16 


12.23 


17 


—88 


16 67 


16 40 


M 2-23 


009 9—09 7 


8 5 


11.58 


11.58 


41 


224 


16 70 


* 


M 2-27 


JJ7.7 Ut. J 


4 8 


12.21 


12.23 




1 7(1 
1 / u 


* 


* 


M 2-31 


006.0-03.6 


5.1 


12.11" 


12.10 


51.0 


157 


* 


* 


M2-33 


002.0-06.2 


5.8 


11.6 a 


11.85 


12.0 


-112 


14 40 


14 40 


M2-39 


008.1-04.7 


3.2 


12.13 


12.07 


8.0 


71 


16 20 


15 80 


M2-42 


008.2-04.8 


3.8 


12.12 


12.10 


14.0 


157 


18.20 


* 


M3-7 


357.1+03.6 


5.8 


12.32 


12.38 


28.0 


-191 


17.30 


16.40 


M3-21 


355.1-06.9 


5.0 d 


11.42 


11.39 


30.0 


-68 


16.20 


15.30 


M3-29 


004.0-11.1 


8.2 


11.7 a 


11.78 


12.0 


50 


15.42 


15.50 


M3-32 


009.4-09.8 


6.0 


11.9° 


11.85 


12.0 


46 


17.40 


17.10 


M 3-33 


009.6-10.6 


5.0 


12.00 


11.93 


7.5 


173 


15.70 


15.90 


NGC 6439 


011.0+05.8 


5.0 


11.71 


11.73 


55.0 


-93 


20.20 


* 


NGC 6565 


003.5-04.6 


13.6 


11.22 


11.25 


38.2 


-20 


* 


18.50 


NGC 6620 


005.8-06.1 


8.0 


11.74 


11.73 


3.5 


72 


* 


19.60 


VY 2-1 


007.0-06.8 


7.0 d 


11.50 


11.57 


37.0 


115 


16.60 




Disk PNe 


















H 1-35 


355.7-03.5 


2.0 


11.52 


11.50 


72.0 


160 


15.70 


15.40 


M 1-29 


359.1-01.7 


7.6 


12.2 a 


12.41 


97.0 


-62 


* 


* 


M 1-61 


019.4-05.3 


1.8 C 


11.43 


11.46 


97.0 


40 


17.10 


* 


NGC 6567 


011.7-00.6 


7.6 


10.95 


10.94 


76.0 


119 


14.42 


14.36 


He 2-118 


327.5+13.3 


5.0 


11.70 


11.67 


10.0 d 


-164 


18.20 


18.70 


IC 4846 


027.6-09.6 


2.0 


11.34 


11.30 


43.0 


151 


15.19 


15.19 



References: A92 - lAckeretalJ ll992|) ; CKS92 - ICatmetalJ Jl992l) ; D98 - 1 Durand et all ll998l) : T89 - iTvlenda et all ll989l) ; T91 - iTvlenda et all Jl99ll) : 

WL07 - from long-slit observations of the current work 

a From lAcker et alN 19921) . 

b Radio measurement from Ziilstra et all )1989|) . 

c Radio measurement from Aaquist & Kwokl ll990l) . 

d Upper limit. 



Table 2. -continued 



Name 


Date 


A-range 


FWHM 


Exp. Time 




(UT) 


(A) 


(A) 


(sec) 


NGC 6567 


07/1995 


3520-7420 


4.5 


20,300 




07/1995 


3900-4980 


0.9 


600,2x1800 




06/2001 


3500-4805 


1.5 


1200 


NGC 6620 


07/1995 


3520-7420 


4.5 


60,300 




07/1995 


4000-4987 


0.9 


2x1800 




06/2001 


3500-4805 


1.5 


3x1800 


VY 2-1 


07/1995 


3520-7420 


4.5 


30,300 




07/1996 


3994-4984 


1.5 


2x1800 



across the nebular surface, thus yielding average spectra for the en- 
tire nebula. 



All the two-dimensional long-slit spectra were reduced us- 
ing the LONG92 package in MIDA^] following the standard pro- 
cedure. Spectra were bias-subtracted, fiat-fielded and cosmic-rays 
removed, and then wavelength calibrated using exposures of a cal- 
ibration lamp. Flux-calibration was carried out using the IRAffl 
package by using wide-slit spectroscopic observations of HST stan- 
dard stars. 

The 3.5 arcmin long-slit was long enough to cover the entire 
nebula for all sample objects while still leaving ample clear area 
to sample the sky background. The sky spectrum was generated 
by choosing sky windows on both sides of the nebular emission 



1 MIDAS is developed and distributed by the European Southern Observa- 
tory. 

2 IRAF is developed and distributed by the National Optical Astronomy 
Observatories. 
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while avoiding regions of background/foreground stars. The sky 
spectrum was then subtracted from the two-dimensional spectrum. 
After sky-subtraction, the two-dimensional spectra were integrated 
along the slit direction over the nebular surface. For a few nebulae 
with a bright central star, such as NGC 6567 and M 2-33, a few 
CCD rows centered on the star were excluded in the summation to 
avoid strong contamination of the stellar light. 

As an example illustrating the high quality of our data, we plot 
in Fig.Q]the spectra of M 3-32 and M 3-21 from 3520 to 4980 A af- 
ter integration along the slit. Also overplotted in Fig.[T]is a synthetic 
spectrum for M 3-21, which includes contribution from recombina- 
tion lines and continua of H I, He I and He II as well as from CELs 
and ORLs emitted by ionized carbon, nitrogen, oxygen and neon 
ions, assuming T c and N e as well as ionic abundances deduced in 
Sections 3 and 4. 

All fluxes were measured on sky-subtracted and extracted one- 
dimensional spectra, using techniques of Gaussian line profile fit- 
ting in MIDAS. For the strongest lines, however, fluxes obtained by 
simply integrating over the observed line profile were adopted. 

H/3 fluxes derived from our wide slit observations are listed 
in Table Q] and compared to values published in the literature. In 
Fig. [2] we compar e our own measurements with those compiled by 
Cahn elai] dl992h . The agreement is quite good except for a few 
nebulae of angular diameters larger than 8". For those large PNe, 
our measurements should be treated as lower limits. 

Most ORLs from CNONe ions are quite weak and often blend 
together (c.f. Fig.[T|l. To retrieve their intensities, multiple Gaussian 
fitting was used. Give n that the ty pical expanding velocity of a PN 
is < 25 km s" 1 (e.g. lKwok|[T994T) , much smaller than the width 
of instrumental broadening (~ 55 km s _1 at H/3 for a FWHM of 
0.9 A, the best resolution of our optical spectra), all lines detected 
with the same instrument setup were thus assumed to have the 
same line width, which was dominated by instrumental broadening. 
This assumption, together with accurately known laboratory wave- 
lengths of identified lines, significantly reduce the number of free 
(non-linear) parameters when fitting blended lines using multiple 
Gaussian. Apart from this procedure, synthetic spectra (c.f. Section 
3.1), were also used to aid line identifications and flux measure- 
ments in cases of serious line blending. The measured line fluxes, 
normalized such that F(H/3) — 100, are presented in Appendix A, 
Table lAll Fluxes of lines measured on high resolution blue spectra 
were normalized to H/3 via H7, the flux of the latter relative to H/3 
was obtained from low resolution spectra which covered nearly the 
whole optical wavelength range. 

Given the weakness of the nebular continuum and high-order 
Balmer lines, and the fact that as n, the principal quantum number 
of the upper level, approaches 20, high-order Balmer lines start to 
blend together, good S/N's and spectral resolution are essential for 
accurate measurements of the Balmer discontinuity and decrement, 
and consequently for accurate determinations of the Balmer jump 
temperature, T C (BJ), and Balmer decrement density, iV c (BD) (c.f. 
Section 3.4). Among our sample of 31 PNe, deep high resolution 
spectra in the blue that cover the wavelength region of the Balmer 
discontinuity and high-order Balmer lines near 3646A are avail- 
able for 8 GBPNe and 3 GDPNe. In subsequent analysis, we shall 
call this sub-sample of PNe as A' and the remaining nebulae (in- 
cluding 17 GBPNe and 3 GDPNe) as sub-sample 'B'. For PNe of 
sub-sample B, values of T C (BJ) and iV (BD) determined are quite 
uncertain. Fortunately, given the fact that emissivities of ORLs have 
only a weak power-law dependence on T c , quite similar to that of 
H I Balmer lines, and are nearly independent on N c at low den- 
sities (N e < 10 6 cm -3 ), the effects of uncertainties in T C (BJ) and 




-12.5 -12 -11.5 -11 

log F(H/S) CKS92 (erg cm" 2 s" 1 ) 



Figure 2. H/3 fluxes derived from our own observations taken with an 8" slit 
width ( ordinate) are comp ared to values published in the literature as com- 
piled by |Cahnetai]h992h (abscissa). The five nebulae for whi ch our mea- 
surem ents are more than 10 per cent lower than those given in lCahn et alj 
1 19921) are marked with open circles and labelled by their names. 



iVc(BD) on ionic abundances deduced from ORLs are minimal (c.f. 
Section 4.2). 

For a few PNe such as M 2-27 in our sample, only low reso- 
lution spectra were available. Given the low spectral resolution and 
S/N ratio, the data were insufficient to deblend and measure weak 
lines. In those cases, synthetic spectra were used to fit the spectra 
in order to obtain crude estimates of T c , N c and abundances from 
ORLs. 

2.2 UV observations 

Sixteen among the 31 PNe have been observed with the IUE in 
the ultraviolet. The data were retrieved from the INES Archive 
Data Server in Vilspa, Spain, processed with the final (NEWSIPS) 
extraction method. All spectra were obtained with the SWP and 
LWP/LWR cameras using the IUE large aperture, a 10.3" x 23" 
oval, large enough to encompass all PNe of the current sample. 
The wavelength coverages of SWP and LWP/LWR spectra are 
from 1150-1975 and from 1910-3300 A, respectively. When sev- 
eral spectra for a given nebula were available, they were co-added 
weighted by integration time. The UV fluxes were normalized to 
F(H/3) = 100 using the absolute H/3 fluxes tabulated in TableQ] 
Absolute H/3 fluxes compiled by CKS92 were adopted, except for 
objects without published reliable measurement, for which our own 
measurements were used (c.f. Section 2.1 and 3.1). A journal of 
IUE observation is given in Table[3] Corrections for interstellar ex- 
tinction will be discussed in Section 3.1. 

Several emission lines from highly-ionized species, such as 
N V A 1240, have been detected in the IUE SWP spectrum of the 
medium ex citation class PN Cn 1 -5. The nebula has a Wolf-Rayet 
central star dAcker & Neineil2003l) . The lines show clearly P-Cygni 
profiles and therefore most likely originate from wind emission 
of the central star. Its optical spectrum yields a He II A4686 in- 
tensity of 0.9 on the scale where F(H/3) = 100. The measure- 
ment was quite uncertain as the 4570 - 4720 A spectral region was 
severely contaminated by strong WR features from C III and C IV. 
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Figure 1. Optical spectra of M 3-32 (upper curve, scaled up by a factor of 2.5) and M 3-21 (middle curve, scaled up by a factor of 1.5) from 3520 - 4980 A, 
showing the rich recombination lines from H I, He I and He II as well as C II, N II, O II, Ne II. Also overplotted is a synthetic spectrum for M 3-21 (lower 
curve), which includes contribution of recombination continuum and line emission from hydrogen and helium as well as CELs and ORLs from CNONe ions. 
The observed spectra have been corrected to laboratory wavelengths using H I Balmer lines and normalized such that 1(11(3) = 100. The spectra have also 
been corrected for interstellar extinction. 
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Table 3. Journal of IUE observations. 



Name 


Data set 


Exp. Time 
(sec) 


Cn 1-5 


SWP 39446 


3000 




SWP 39455 


6000 




LWP 18567 


4000 




LWP 18577 


3000 


H 1-42 


SWP 39447 


12840 


H 1-50 


SWP 54686 


23100 


IC 4846 


SWP 33381 


3600 




LWP 13130 


4800 


M2-23 


SWP 44207 


25380 


M 2-33 


SWP 55509 


24000 


ivi z-jy 




ojoU 


M3-21 


SWP 39464 


8460 




SWP 39465 


10800 


M3-29 


SWP 30480 


1500 




SWP 54464 


8400 


M3-32 


SWP 54726 


24000 




SWP 55541 


24000 


M 3-33 


SWP 44199 


7200 


JNOL. 04jy 


SWr Jj454 


t yzuu 


NGC 6565 


SWP 07984 


2100 




SWP 35690 


9000 




SWP 24266 


3600 




LWr 040 i 1 


3600 




LWR 06955 


2100 




LWP 15144 


9000 


NGC 6567 


SWP 17019 


11160 




SWP 45362 


7200 




SWP 45381 


1320 




LWR 13307 


3600 




LWP 23708 


3600 


NGC 6620 


SWP 45836 


1800 




SWP 48320 


8760 




SWP 54512 


2400 




SWP 54700 


3300 


Vy 2-1 


SWP 44200 


12240 



At T c = 10 000 K, the He II A1640/A4686 intensity ratio has a pre- 
dicted value of 6.6. The measured optical flux of the He II A4686 
line then implies a A1640 line flux of 6, which is a factor of ten 
lower than measured from the IUE spectrum (after reddening cor- 
rections). This shows that He II emission detected in the IUE spec- 
trum is indeed dominated by wind emission, rather than from the 
gaseous envelope. In Section 4, we will use the observed intensity 
of the He II A4686 to obtain an estimate of the He ++ /H + abun- 
dance for the gaseous envelope of this PN. We also note that the 
blue absorption is quite weak for the C IV A1550 feature and its ob- 
served flux will be used to obtain a rough estimate of the C 3+ /H + 
ionic abundance for this PN. 

2.3 Infrared observations 

Three objects in our sample were observed with the Short Wave- 
length Spectrometer (SWS) and Long Wavelength Spectrometer 
(LWS) on board ISO. Cn 1-5 has one SWS01 and one LWS01 spec- 
trum available. The SWS01 spectrum covered the 2.4-45 /im wave- 



length range at a FWHM of approximately 0.3 /im, whereas the 
LWS01 covered from 43-197 fxm at a FWHM of about 0.6 /im. For 
M 2-23, several lines in the SWS wavelength range were scanned 
using the SWS02 mode, which yielded a spectral resolving power 
between 1000 and 4000, depending on wavelength as well as on 
source angular extension. Finally the SWS02 and LWS02 observ- 
ing modes were used to scan a number of lines in NGC 6567 in 
both the SWS and LWS wavelength ranges. 

The angular sizes of the three nebulae are small enough such 
that both SWS and LWS observations should capture the total flux 
from the entire nebula. Line fluxes measured from the infrared 
spectra were normalized to F(H(3) = 100 via the absolute H/3 flux 
listed in Table[T]which were compiled by CKS92. The normalized 
line fluxes are listed in Table lA2l and then dereddened as discussed 
in the next Section. 



3 NEBULAR ANALYSIS 

Our a nalyses follow closely the procedures outlined in iLiu et all 
< l200ll) who presented detailed studies of two GBPNe, M 1-42 
and M 2-36. Firstly, extinction curves towards individual sightlines 
were derived from radio continuum flux density and H I, He I and 
He II recombination lines. Plasma diagnostic analyses were then 
carried out using both CELs and ORLs, followed by ionic and total 
elemental abundance determinations. 

3.1 Extinction correction 

Before any nebular analyses, measured line fluxes need to be cor- 
rected for effects of extinction by intervening dust grains along the 
sightline. Some PNe, such as the nearby NGC 7027 and NG C 6302 
are a l so known to harbour large amounts of loca l dust (Seaton 
1 19791 : iMiddlemassll 19901 : iLester & Dinersteirjl 19841) . For GDPNe, 
it is generally sufficient to use the standard Galactic extinction law 
for the general diffuse ISM, which has a total-to-selective extinc- 
tion ratio_Rv = A(V)/75 B -v = A(V)/[A(B) - A(V)] = 3.1 
(c.f. lHowartr]| 19831 : H83 hereafter). In that case, the amount of ex- 
tinction can be characterized by a single parameter, the logarithmic 
extinction at H/3, c(H/3). However, large variations in reddening 
curves tow ards individual sightlin es are also observed, particularly 
in the UV |m assa & Savage! 19891) . It is found that the linear rise of 
extinction in the far-UV decreases dramatically with R v and sug- 
gests that in dense regions, carriers of the far-UV exti nction, includ- 
ing those of the 2175 A bump reside in large grains ( Carde lli et al.l 
Il989l) . Earlier work on extinction towards the Galactic bulge in the 
optical and UV region indicates a steep redde ning law in the UV , 
with values of R v ranging from 1.75 to 2.7 JWalton et alJl 19931 : 
iRuffle et al.ll2004h. 

Liu et al. I feOOll) derived extinction curves towards M 1-42 
and M 2-36 by comparing the observed Balmer decrement, the 
He II A1640/A4686 and A3203/A4686 ratios, and the ratio of to- 
tal H/3 flux to radio f-f co ntinuum flux density wi th the predictions 
of recombination theory l lStorev & Humme Il995l) . In the current 
work the same method was adopted. Firstly, intrinsic intensities 
7(A) of H I Balmer lines were calculated from the radio contin- 
uum flux density at 5 GHz or 1 .4 GHz, assuming that the nebula is 
optically thin at those frequencies. Comparison between the pre- 
dicted and observed fluxes of individual lines yields the extinction 
value A(X) at the wavelength of the line. Similarly, from the ob- 
served He IIA1640/4686 ratio, the differential extinction between 
the two wavelengths, y4(A1640) — A(A4686), can be derived. The 
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value, when combined with yl(A4686), obtained by interpolating 
A(X4340) and vl(A4861), then yields A(A1640). In this way, ex- 
tinction values at a set of wavelengths can be established, allowing 
characterization of the extinction law towards individual objects. 

Since the ratio of the radio continuum flux density to F(H/3) 
depends weakly on T c and helium ionic abundances, a prelimi- 
nary plasma diagnostic and abundance anal ysis were carrie d out 
assuming the standard ISM reddening law l lHowartf]|l983l) . Ra- 
dio flux de nsities at 1.4 GHz for ou r sample nebulae were mostly 
taken from lCondon & Kaplan! Tl 998). whereas those of 5 GHz were 
mainly from Iziilstra et alj 1 1989l).For total H /3 flux, F(B.0), we 
adopted values compiled by Cah n et alj dl992h . For objects without 
published reliable measurement of F(H/3), we have used our own 
measurements obtained using an 8-arcsec wide slit. We check the 
accuracy of our measurements in Fig.[2]by comparing our wide-slit 
measurements and those from the literature. The agreement is good. 
The four out of five objects for which our measured values are over 
10 per cent smaller than those given in the literature have optical di- 
ameters larger than 7" and thus the discrepancies are likely caused 
by our finite slit width. Fo r M 2-33, a compa ct nebula with an opti- 
cal diameter less than 6" (Dopita et al. 199 (j), we prefer to us e our 
own measurement as the value given by lAcker et al.1 J 19921) was 
noted to be quite uncertain. 

From the radio flux and total H/3 flux and adopting preliminary 
values of derived T c and He ionic abundances, we obtain A(H/3), 
the absolute extinction at H/3, which was then used together with 
the observed intensity ratios of H I, He I and He II lines to calcu- 
late extinction values at the corresponding wavelength. During this 
process, some cautions have to be exercised. For example, in some 
PNe, He was blueshifted to the position of the interstellar Ca II H 
line and absorbed by it. Similarly, in some PNe the He II A 1640 
line is contaminated by wind emission from the central star. 

As an example, Fig.[3]plots X towards NGC 6565 as a func- 
tion of x, where x = l/\(pm) and X(x) = E(X - V)/E(B - 
V) = [.4(A) - A(V)]/A{B) - A(V)]. Also overplotted i n the Fig- 
ure ar e standard Galactic IS M extinction laws gi ven by [Howartb 
dl983l) (H83; dotted line) and lCardelli e"tai]dl989h (CCM89: long- 
dashed line) for a total-to-selective ratio Ry = 3.1. The extinc- 
tion curve deduced for the sightline towards NGC 656 5 is quite 
simila r to previously found for the bulge PN M 1-42 by iLiu et all 
j200lh - the extinction in the far ultraviolet as derived from the 
He II A1640/A4686 ratio is much higher than predicted by the 
standard Galactic ISM extinction law. It is unfortunate that there 
are no lines usable near the promin ent 2175 A bum p present in the 
standard extinction law. Following ILiu et aT 1 42001b, we fit X as a 
function of x with a straight line, which seems to be appropriate for 
available data. For NGC 6565, the result is also shown in Fig. [3] 

Excluding He 2-118, for which only a lower limit on its 
5 GHz flux density is available from the literature, for the remain- 
ing 30 sample PNe, the normalized extinction, X(x), was derived 
for each object at wavelengths of all usable lines. The logarith- 
mic extinction, /(A), commonly used in nebular analyses and de- 
fined such that /(A) = at H/3, was also calculated. /(A) = 
0.4A(A)/c(H/3) - 1, where c(H/3) = 0.4A(H/3) is the logarith- 
mic extinction at H/3. The deduced normalized extinction curve was 
then compared with the H83 law and the CCM89 law. For eight 
sample nebulae, reliable fluxes of the He II A1640 line and/or the 
A3203 line are available, providing an estimate of the extinction in 
the UV. For four of them, the extinction X(x) shows a very good 
linear relation extending from optical to the UV. For the four ob- 
jects, a linear function was used to fit X(x) as a function of x. The 
resultant extinction curve was then used to deredden the observed 
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Figure 3. Normalized extinction X(x) towards NGC 6565. The solid line 
is a linear fit to the data, X(x) = —2.69 + 1.34a;. The dotted line and 
long-dashed lines are standard Galactic reddening law for Ry =3.1 given 
by H83 and CCM89, respectively. 

line flux. For the other four PNe, either the H83 or the CCM89 law 
was used to fit the data, depending on which one yielded a better 
fit. We note that the two extinction laws have noticeable differences 
in the 3500-4500 A wavelength range. 

For the 22 remaining nebulae for which no estimate of UV ex- 
tinction was possible, the extinction curve was either fitted with a 
linear function, the H83 law or the CCM89 law, whichever yielded 
the best fit. The extinction law used for each nebula is listed in the 
last column of Table [4] If an entry in this column consists of two 
numbers, it denotes a linear extinction law, fix) = a + bx, was 
used, and the two numbers refer to fitted values of parameters a 
and b, respectively. In Table|4] we also list values of c(H/3) derived 
from ratios of Ha/H/3 (Column 2), F„(5GHz)/H/3 (Column 3) and 
of He II A1640/A4686 (Column 4), as well as our final adopted val- 
ues (Column 5) used to deredden the observed line fluxes. Finally, 
Column6ofTableH gives the total-to-selective extinction ratio Ry. 



3.2 CEL plasma diagnostics 

Accurate determinations of abundances rely on the reliable mea- 
surements of T e and N e . In what follows, we present de- 
tailed plasma diagnostics using both CELs and recombination 
lines/continua. Values of N c and T c , have been derived from CEL 
diagnostic ratios by solving equations of statistical equilibrium of 
multi-level (^ 5) atomic models using EQUIB, a Fortran code orig- 
inally written by I. Howarth and S. Adams. CEL diagnostic lines 
used are tabulated in Table [5] together with the ionization poten- 
tials required to create the emitting ions. The atomic p arameters 
used i n the present work are the same as those used by ILiu et al.1 
l l2000h in their case study of NGC 6153. 

In order to obtain self-consistent results for N D and T c probed 
by CEL lines, we proceed as following. Firstly, iV c ([Cl III]) and/or 
N c ([Ai IV]) were obtained assuming a T c of 10000 K and then 
the average value of the derived iV c was adopted in calculating 
T c ([0 III]). The calculations were iterated to obtain final self- 
consistent values of N e and T c . We assume the results thus ob- 
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Table 6. Electron temperatures and densities from CELs. 'H' and 'L' denote the measured line ratio exceeds its high and low temperature/density limit, 
respectively, whereas 'E' indicates that recombination contribution as estimated using Eqs. (1) and (2) exceeds the actual measured flux of the [N II] A5755 
line or of the [O II] AA7320, 7330 lines. 



Nebula 


[0 m] na 


T c (CELs) (K) 
[N ll] na 


[O n] na 


He 2-118 


12630 


14960 


H 


M2-4 


8570 


9920 


14570 


M2-6 


10100 


9730 


11650 


M3-7 


7670 


6900 


5110 


M 1-20 


9860 


11180 


16570 


NGC 6439 


10360 


9270 


11340 


H 1-35 


9060 


12080 


13720 


M 1-29 


10830 


9020 


9000 


Cn 2-1 


10250 


12030 


H 


H 1-41 


9800 


9530 


12790 


H 1-42 


9690 


9050 


13270 


M2-23 


11980 


H 


H 


M3-21 


9790 


12800 


19000 


H 1-50 


10950 


12070 


15360 


M2-27 


11980 


8650 


8840 


H 1-54 


9540 


11600 


11360 


NGC 6565 


10300 


10100 


10200 


M2-31 


9840 


11370 


13470 


NGC 6567 


10580 


10016 


14360 


M2-33 


8040 


9150 


11070 


IC 4699 


11720 


12490 


9390 


M2-39 


8050 


8300 


11030 


M2-42 


8470 


9350 


11860 


NGC 6620 


9590 


8630 


9630 


Vy 2-1 


7860 


8580 


8750 


Cn 1-5 


8770 


8250 


10330 


M3-29 


9190 


8750 


8930 


M3-32 


8860 


8230 


E 


M 1-61 


8900 


11510 


8460 


M3-33 


10380 




7480 


IC 4846 


9930 


13280 


H 



AT e (CELs) (cm" 3 ) 

[OlI]nn [Sll]nn [CI III]nn [Ar IV] 
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Figure 4. CEL plasma diagnostic diagram for NGC 6567. For diagnostic 
ID's, c.f. Tabled 



Figure 5. Plasma diagnostic diagram for M 2-23. For diagnostic ID's, c.f. 
Table \5\ The four broken lines represent density diagnostics from [Fe III] 
lines (c.f. text for more details). 
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Table 4. Extinction towards individual nebulae. 



















opt. 


rad. 


UV 


adopt 






Cn 1-5 


0.49 


0.39 


0.59 b 


0.49 


3.05 


-0.91,0.44 


Cn 2-1 


1.07 


0.83 




0.83 


3.10 


CCM89 


H 1-35 


1.51 


0.96 




1.31 


3.27 


CCM89 


H 1-41 


0.65 


0.67 




0.65 


3.45 


-1.21,0.59 


H 1-42 


0.87 


0.74 




0.77 


3.35 


-1.31,0.64 


H 1-50 


0.68 


0.53 




0.68 


3.70 


-1.15,0.56 


H 1-54 


1.54 


0.68 




1.53 


1.92 


-2.26,1.09 


He 2-118" 


0.17 


0.14 




0.16 


3.85 




IC 4699 


0.23 


0.45 


0.16 


0.22 


3.21 


-1.42,0.69 


IC 4846 


0.69 


0.34 




0.65 


2.79 


CCM89 


M 1-20 


1.40 


1.02 




1.13 


3.40 


CCM89 


M 1-29 


2.03 


1.63 




1.83 


3.10 


CCM89 


M 1-61 


0.92 


0.92 




0.92 


1.81 


-2.31,1.12 


M2-4 


1.33 


0.81 




0.81 


1.97 


CCM89 


M2-6 


1.14 


0.84 




0.85 


2.62 


CCM89 


M2-23 


1.20 


0.98 




1.12 


2.97 


CCM89 


M2-27 


1.31 


1.37 




1.35 


2.61 


CCM89 


M2-31 


1.41 


1.18 




1.15 


3.07 


CCM89 


M2-33 


0.55 


0.61 




0.55 


3.00 


CCM89 


M2-39 


0.61 


0.63 




0.63 


1.92 


CCM89 


M2-42 


1.06 


0.66 




0.65 


1.97 


CCM89 


M3-7 


1.65 


1.51 




1.46 


4.18 


H83 


M3-21 


0.50 


0.37 


0.78 


0.35 


2.92 


-1.27,0.62 


M3-29 


0.24 


0.12 




0.12 


1.86 


CCM89 


M3-32 


0.64 


0.40 


0.69 


0.62 


3.10 


H83 


M3-33 


0.50 


0.42 


0.32 


0.39 


3.31 


CCM89 


N6439 


1.10 


0.93 


1.15 


0.93 


2.81 


CCM89 


N6565 


0.32 


0.36 


0.52 


0.37 


3.10 


-1.03,0.50 


N6567 


0.90 


0.66 


().73 c 


0.73 


3.36 


H83 


N6620 


0.52 


0.46 


0.86 


0.43 


3.10 


CCM89 


Vy 2-1 


0.83 


0.66 




0.64 


2.80 


CCM89 



a Only a lower limit for the 5 GHz flux is available from the literature. 
Thus c(H/3) derived from the radio continuum to H/3 flux ratio is an upper 
limit, as is Ry ■ 

6 Derived from the He II A1640/A3203 ratio. 
c Derived from the UV continuum. 

d Extinction law used to fit the extinction curve. H83 and CCM 89 de- 
note extinction laws given by Howarth ( 1983) and Cardelli et al. ( 1989), 
respectively. When two numbers are listed, it denotes a linear extinction 
law, /(A) = a + bx, and the two numbers refer to values of a and b, 
respectively, x = l/A(/im). 

tained represent physical conditions of the high ionization regions. 
Similarly, T ([N II]) was derived in combination with iV c ([0 II]) 
and 7Y C ([S II]) to represent the physical conditions in the low ion- 
ization zones. If all observed optical 7V C -indicators yielded similar 
values of N c , the results were averaged and used for the whole neb- 
ulae, both the high and low ionized regions. Recombination contri- 
butions to intensities of the [O II] and [N II] auroral lines were 
corrected for using the formulae derived by iLiu et al] d2000h . 



/r(A5755) _ olm0 .3o v N 2 + 
____ _ 6 . wt x — ; 

7 R (A7320 + A7330) 



n ocj-0.44 O 
= 9.36t X 



2+ 

H~ ! 



(1) 



(2) 



where 2+ /H+ and N^/H" 1 
ORLs, and t ee T C /10 4 K. 



Table 5. Plasma diagnostic line ratios and ionization potentials of the emit- 
ting ions. Subscripts 'n', 'a' and *f represent nebular, auroral and fine- 
structure lines, respectively. 



in 
ID 




lr (ev) 




T c -diagnostics 




[OlI] na 


/(A7320 + A7330)//(A3726 + A3729) 


13.6 


[NH] na 


/(A6548 + A6584)//(A5754) 


14.5 


[O IIl]na 


/(A4959 + A5007)//(A4363) 


35.1 


[Nem] nf 


7(15.5/mi)//(A3868) 


41.0 


[ArIII] nf 


7(9.0iim)/7(A7135) 


27.6 




-/V e -diagnostics 




[Sll] nn 


7(A6731)//(A6716) 


10.4 


[On] nn 


7(A3729)//(A3726) 


13.6 


[CllII]„„ 


7(A5537)//(A5517) 


23.8 


[ArIV] nn 


7(A4740)/7(A4711) 


40.7 


[Nem] ff 


7(15.5Mm)//(36.0^m) 


41.0 



IWang et al] te004h have examined the four optical 7V e - 
diagnostics, [O II] A3729/A3726, [S II] A6716/A6731, 
[CI III] A5517/A5537 and [Ar IV] A4711/A4740, for a sam- 
ple of over 100 PNe, including 31 PNe in the current sample. The 
atomic parameters recommended in their work were adopted in the 
current study. For objects in common, N c presented here differ 
by no more than 1 per cent from those published by IWang et all 
d2004l) . 

Fig. [4] and [5] present plasma diagnostic diagrams for 
NGC 6567 and M 2-23, respectively, as two examples. In the case 
of NGC 6567, various 7V C -diagnostic lines yield consistent values 
of 7Y e , so a constant N c has been adopted in calculating T c from 
the [N Il] na and [O IIl] n a nebular to auroral diagnostic line ratios. 
The two To's thus obtained were then adopted in calculating CEL 
abundances for singly ionized species and for species of higher ion- 
ization degrees, respectively. 

Electron temperatures and densities derived from various CEL 
diagnostics are tabulated in Table [6] In the Table, 'H' and 'L' de- 
note the measured line ratio exceeds its high and low tempera- 
ture/density limit, respectively, whereas 'E' indicates that recom- 
bination contribution as estimated using Eqs.(l) and (2) exceeds 
the actual measured flux of the [N II] A5755 line or of the [O II] 
AA7320, 7330 lines. 

Kingsburgh & Barlow (1994, KB94 hereafter) have examined 
the relation between T c ([0 III]) and T C ([N II]) for a sample of 
Galactic PNe. They found that the ratio of the former to the later 
increases with increasing /(A4686), the intensity of He II A4686, 



T C ([NII]) 



= (1.15 ± 0.06) + (0.0037 ± 0.0009)7(A4686), 



with a correlation coefficient of 0.69. A similar analysis for GBPNe 
in the WL07 sample and GDPNe in the TLW sample (c.f. §[7} ex- 
hibits similar positive correlation, as shown in Fig. [6] A linear fit to 
the 45 points in the Figure yields: 



r c ([Oin]) 

T C ([NII]) 



(0.91 ± 0.02) + (0.0049 ± 0.0008)J(A4686), 



are ionic abundances derived from 



with a correlation coe fficient of 0.70, which i s roug hly in agreement 
with that obtained by Kingsburg h & Barlowl dl994l) . 
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Table 7. Electron temperatures deduced from H I Balmer discontinuity, from He I recombination lines and from O II ORLs, as well as electron density derived 
from the H I Balmer decrement. An 'L' denotes that the measured O II A4089/A4649 exceeds its low temperature limit of 0.41 at 288 K. 



Nebula T C (H I) (K) T c (He I) (K) T c (OlI)(K) 7V C (H I) (cm -3 ) 

RT A 66 78 A5876 A7281 A7281 A4649 pjp, 

J A4471 A4471 A5876 A6678 A4089 



He 2-118 


14500 


12050 


12870 


5250 


6800 




5.0 


M2-4 


7900 




14100 


6100 


6100 




4.6 


M2-6 


11700 


4300 


2860 


6200 


6500 


L 




M3-7 


6900 


8100 


2580 


1600 


2500 






M 1-20 


12000 






5100 


6000 




4.2 


NGC 6439 


9900 


2060 


1810 


4810 


4900 


851 


5.5 


H 1-35 


12000 






5800 


7300 


334 


4.8 


M 1-29 


10000 


6900 




4000 


4100 






Cn 2-1 


10800 


6750 




5050 


5450 


857 


5.0 


H 1-41 


4500 


7110 


5770 


2690 


2930 


L 




H 1-42 


10000 


3700 


2250 


5890 


6290 


662 




M 2-23 


5000 


4900 


1320 


6100 


7330 


410 




M 3-21 


10400 


3480 


2150 


5170 


5520 


602 


5.0 


H 1-50 


12000 


9050 


3100 


5420 


6480 


L 


5.1 


M2-27 


14000 


8410 




2280 


2920 






H 1-54 


12500 


13600 


10400 


4270 


5840 






NGC 6565 


8500 


3850 


6400 


4620 


4090 


11722 




M2-31 


14000 


2170 


1220 


3820 


4460 






NGC 6567 


14000 


7850 


2670 


6720 


7790 




4.0 


M2-33 


7000 


8510 


5200 


4050 


4650 


L 




IC 4699 


12000 


7300 


5950 


3310 


2460 


L 




M2-39 


5500 


6810 


7070 


6610 


6600 


L 




M2-42 


14000 


8000 


7700 


3600 


3300 






NGC 6620 


8200 


2500 


2400 


3800 


3800 


3162 


4.0 


Vy 2-1 


8700 


5500 


4550 


4750 


4970 


L 




Cn 1-5 


10000 


11600 




2450 


2940 


18000 




M3-29 


10700 


6800 


4750 


1600 


1800 






M 3-32 


4400 


3860 


3100 


1560 


1710 


L 


3.6 


M 1-61 


9500 


18500 




3500 


5490 






M3-33 


5900 


6650 


3940 


4380 


5020 


1465 




IC 4846 


18000 


6530 


2430 


6080 


6970 


9954 





3.3 Two interesting objects 

For most nebulae in our sample, diagnostics that sample respec- 
tively the high and the low ionized regions generally yield N e of 
relatively small differences, as in the case of NGC 6567 (Fig.gJ. 
There are however a few exceptions, such as M 2-23. For this ob- 
ject, as shown in Fig. [5] lines emitted by ions of high ionization po- 
tentials (i.e., [Ar iv] nn and [Ne m]ff) yield a A^ of ~ 10 4 7 cm -3 , 
approximately 0.8 dex higher than derived from lines of low ion- 
ization potentials, such as [S Il] Iln . We note that for this nebula, the 
[N Il] na and [O Il] lm nebular to auroral line ratios become sensitive 
to N c as well as to T e . M 2-23 thus shares some similarities with 
Mz 3 and M 2-24, in that both have been previously found to har- 
bour a dense emissio n core surrounded by an o uter lobe of much 
lower densities (c.f. IZhang & LiJ 120021 l2003h . In IZhang & LiJ 
(2002), high critical density [Fe III] diagnostic lines were used to 
probe the physical conditions in the dense central emission core of 
Mz 3. The [Fe III] line strengths observed in M 2-23 are on the av- 
erage about a third of those detected in Mz 3. Unfortunately, lines 
sensitive to T e , such as the A6096 or the A7088 lines were all too 
faint to be detectable in our current spectra of M 2-23. In Fig. [5] 
we show, from left to right, the loci of four [Fe III] Ac-sensitive di- 
agnostic line ratios, A4659/A4770, A4755/A4881, A4770/A4778 
and A4702/A4734. Except for the first line ratio, which is less reli- 
able given the intrinsic weakness of the [Fe III] A4770 line, all the 



other three consistently yield a N e ~ 10 6 cm 3 , comparable to the 
value of ~ 10 6 - 5 cm" 3 found for Mz 3 bv lZhang & LiJ d2002[) and 
considerably higher than values given by other N c -diagnostics of 
lower critical densities. 

In the ISO SWS spectrum of M 2-23, apart from the ionic 
fine-structure lines from heavy elements, H I Bra, Br/3 and Pfa 
lines were also well detected. The relative strengths of those high 
excitation H I recombination lines have a weak dependent on T e . 
Using the emission coefficients given bv lStorev & Hummer! fl995h . 
we determined a T a of 5 000 K, with an uncertainty about 1 500 K. 
The value is very similar to that deduced from the H I Balmer dis- 
continuity. 

Another object in our sample that exhibits some interesting 
characteristics is M 2-39. For this nebula, all Ac -diagnostic line ra- 
tios measurable in the optical, [S Il] n n, [CI IIl]nn and [Ar IV] nn , 
consistently yield a low N c between 1500 and 5100 cm -3 , 
whereas the [N Il] na ratio gives a T c of 8 300 K (Table |6). Yet 
the observed [O IIl] n a ratio shows an abnormally low value of 25, 
indicating a T c of about 28 000 K for a N e of below 10 000 cm -3 . 
The [O IIl] na temperature could be lower if the density prevailing 
in the 2+ zone is much higher than 10 000cm -3 . Some evidence 
in favour of a dense 2+ zone is provided by a well observed fea- 
ture near 7171 A, which we tentatively identified as the [Ar IV] 
A7 170.62 line. If the identification is correct, then its intensity rela- 
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Figure 6. T e ([0 III])/T C ([N II]) ratios against /(A4686) for GBPNe (filled 
circles) and GDPNe (open circles). The solid line is a linear fit to the 45 
points in the Figure. 



tive to the [Ar IV] A A471 1,4740 nebular lines would indicate a N e 
of 10 6 ' 6 cm -3 . For comparison, if we assume a T c of 10000 K, as 
indicated by [N Il] n a and [O Il]na, then a iV c as high as 10 6 ' 2 cm" 3 
would be required to reproduce the measured [O III] na ratio. The 
identification of the feature at 7171 A as [Ar IV] was however 
marred by the fact that another [Ar IV] line, at 7263 A, expected 
to have a strength only 30 per cent lower than that of the [Ar IV] 
A7171, was either absent or had at most a strength only one third 
of the measured strength of the 7171 A feature. Clearly, better data 
are needed to clarify whether M2-39 also contains a dense high 
excitation emission core as in the case of Mz 3, M 2-24 and proba- 
bly also of M 2-23. In our current analysis, we have adopted a N c 
of 10 6 ' 6 cm -3 for the high ionization regions of M 2-39, assuming 
that the 7171 A is indeed entirely due to [Ar IV]. Then the observed 
[O IIl]na ratio yields a T c of 8 050 K (c.f. Table©. This N c and 
T e have been used to calculate ionic abundances of heavy element 
ions from the high ionization regions. We note that the abundances 
thus obtained were abnormally high compared to other bugle PNe 
in the sample, even higher than values deduced from ORLs for the 
same object in the case of carbon and oxygen (c.f. Tablell4t. These 
abundances are however quite uncertain and should be treated with 
caution. 

Given their peculiarities, both M 2-23 and M 2-39 will be 
excluded in our statistical analysis of the abundance patterns of 
GBPNe, leaving a sample size of 23. 



3.4 Hydrogen recombination temperatures and densities 

Electron temperatures have been determined for all sample PNe 
from the Balmer discontinuity of H I recombination spectra with 
various levels of uncertainty depending on the resolution and S/N 
ratios of spectra available for the wavelength range of the Balmer 
ju mp at 3650 A. Fo r practical reasons of measurement and follow- 
ing lLiuetal.l f2001), we define the Balmer jump, BJ, as the differ- 
ence of continuum flux densities measured at 3643 and 3681 A, i.e. 
J C (A3643) — J C (A3681). Balmer jump temperature, T C (BJ), was 
then calculated from the ratio of Ba lmer jump to H 1 1 using the 
formula provided bv lliu et"aT] ( feOOll) . 




6 X [I(X7281)/[(X6678)] 



- 



1000 10000 
T e (K) 

Figure 7. He I line ratios as a function of T c , assuming N c = 4000cm -3 . 
The ratios observed in M 3-32 are marked. 

R 1 3 

T c = 368 x (l + 0.259y+ +3.409y 2+ )( — (3) 

where BJ/H11 is in units of A" 1 and y + and y 2+ denote he- 
lium ionic abundances, He + /H + and He 2+ /H + , respectively. Since 
y + and y 2+ derived from helium recombination lines also have a 
weak dependence on adopted T c , the process was iterated until self- 
consistent values for y + and y 2+ and T C (BJ) were achieved. The 
results are tabulated in the second column of Table|7] 

Electron densities, 7Vc(BD)'s, were also derived by comparing 
the observed decrement of high-order Balmer lines with that pre- 
dicted by recombination theory and are presented in the last column 
ofT able|7l 

IZhang et alj d2004l) present determinations of T C (BJ) and 
N C (BD) for a sample of 48 Galactic PNe, including some objects 
in the current sample, by fitting the observed spectra with synthe- 
sized theoretical spectra with T a and N c as input parameters. The 
technique only works well for high resolution spectra. For objects 
in common, our current results are in good agreement with those 
presented in Zh ang et alj ([2004). For objects having high resolu- 
tion spectra, typical uncertainties in T C (BJ) are about 1000 K. For 
those for which only low resolution spectra are available, the un- 
certainties increase to about 2000 K. 

3.5 Helium temperatures 

Electron temperatures have been derived from He I recombi- 
nation line ratios A6678/A4471, A5876/A4471, A7281/A5876, 
A7281/A6678 and are presented in Cols. 3-6 of Table [7] Albeit 
the lines involved are weaker, T a determined from the fourth pair 
of diagnostic line ratio has a number of advantages ove r those de- 
duced from the other three pairs (c.f. Zhang et al. 2005 for details) 
and will be adopted in the current worio In Fig. [7] we show vari- 
ations of the four He I line ratios as a function of T a for an as- 
sumed N c = 4000 cm -3 . Values measured in M 3-32 are marked as 
an example. The parameters us ed to calculate em issivities of He I 
lin es are taken f r om Ta ble 2 of IZhang et al.1 d2005l) . In their analy- 
sis, |zhang|etal] (2005) presented He I temperatures deduced from 

3 However, We note that the A728 1 line is quite sensitive to the assumption 
of Case A or B recombination. Although Case B is generally regarded as a 
good approximation, there is some evidence of departure from pure Case B 
recombination (Liu et al. 2000, 2001). In that case, T c (HeI) deduced from 
the A7281/A6678 and A7281/A5876 ratio would be underestimated 
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the A7281/A5876 and A7281/A6678 ratios for a large sample of 
Galactic PNe; 9 of them are also in the current sample. For ob- 
jects in common, T c obtained by them are slightly higher than ours. 
The differences are mainly caused by the differences in extinction 
corrections. In their analysis, they had adopted the standard H83 
extinction law for the general ISM and determined extinction to in- 
dividual nebulae based on optical data alone. In the current work, 
a much more detailed and thorough treatment of the extinction has 
been carried out (c.f. Section 3.1). 

3.6 Temperatures from heavy element ORLs 

iLiu et alj ( l200lh compared the observed relative intensities of O II 
ORLs with predictions of the recombination theory in four PNe ex- 
hibiting particularly large adfs, M 1-42, M 2-36, NGC 7009 and 
NGC 6153. They found that while the agreement was quite good 
in general, there were exceptions. For example, the intensity of the 
strongest 3p-3s transition, 3p 4 D^ , 2 - 3s 4 P 5 / 2 A4649, appears to be 
too strong (by ~ 40 per cent) relative to the strongest 4f - 3d tran- 
sition, 4fG[5]°jy 2 - 3d 4 F 9 / 2 A4089, compared to the theoretical 
value calculated at T c (BJ). It was later realized that the discrep- 
ancies were actually caused by the fact that those heavy element 
ORLs arise mainly from another component of plasma of very low 
T e , much lower than values dedu ced from the Balmer discontinuity 
of H I recombination spectrum (Liu 2003). By comparing observed 
ratios of the two strongest O II ORLs, A4089/A4649, with the the- 
oretical predictions calculated down to a minimum T c of 288 K, 
values of T c (0 II), the average T c under which O II OR Ls are emit- 
ted, have been determined for a large samp le of PNe dTsamis et al.1 
l2004lLiu etai]|2004bl.lWesson et alj2005l) . 

For the current sample of PNe, we have been able to determine 
T e (0 II) for 1 1 objects from the A4089/A4649 ratio. The results are 
presented in Col. 7 of Table|7] For another 8 nebulae in the sample 
(flagged with an 'L' in Tabled for which both lines have been de- 
tected, the measured line ratios exceed 0.41, the maximum value 
at the low T c limit of 288 K. Part of the discrepancy was clearly 
caused by measurement uncertainties, given the weakness of the 
lines. There is also an opportunity that the O II A4089.3 may have 
been contaminated by the Si IV A4088.8 line, lead ing to artificially 
high A4089/A4649 ratios in some PNe jLiull2006bl: c.f. also Secc- 
tion6.1). 



4 IONIC AND TOTAL ELEMENTAL ABUNDANCES 
4.1 Ionic abundances from CELs 

Abundances for ionic species of C, N, O, Ne, S, Ar, CI and Fe have 
been determined from CELs for 31 sample nebulae. In the calcula- 
tions, we have adopted T a derived from [N Il] na nebular to auroral 
line ratio and N c derived from the [S II] nn nebular line ratio (or the 
average of values derived from the [S Il] nn and [O Il] nn ratios if 
the latter is available). For ionic species of higher ionization stages, 
T e deduced from the [O IIl] na ratio and the average iV c derived 
from the [CI IIl] nn and [Ar IV] nn ratios were used. Note that for 
most objects, we have corrected for a recombination contribution 
to the [N II] A5755 auroral line, except for a few cases where the 
contamination as estimated using the aforementioned formula (and 
adopting an ionic abundance N 2+ /H + as yielded by N II ORLs) 
was found to be unphysically large, becoming comparable to or 
even larger than the actually measured total flux of the A5755 fea- 
ture. In the case of M 3-33, the [N II] A5755 line was too faint to 



have a reliable measurement. For this nebula T c [0 III] has been 
used to calculate abundances of all ionic species. 

Ionic abundance ratios, X 1+ /H + , derived from collisionally 
excited UV, optical and IR lines are tabulated in Tables[8](oxygen), 
l9l (carbon). 1101 (nitrogen). [TT1 (neon and argon) and 1 121 (sulfur and 
chlorine). 

4.2 Ionic abundances from ORLs 

4.2.1 He 

Singly and doubly ionized helium abundances were derived from 
the He I A4471, A5876 and A6678 recombination lines and from 
the He II A4686 recombination line, respectively. The results from 
the three He I lines were averaged, weighted by 1:3:1. The total 
helium abundances relative to hydrogen were obtained by summing 
abundances of singly and doubly ionized species. The results are 
tabulated in Cols. 2-4 of Tabled 

4.2.2 C, N, O, Ne and Mg 

Ionic and elemental abundances have been derived from ORLs for 
some or all of the elements C, N, O, Ne and Mg for our sample of 
31 PNe. Electron temperatures derived from the H I Balmer dis- 
continuity, T C (BJ), were adopted in the calculations. For the 11 
PNe of sub-sample A, A r c (BD) were derivable and were adopted. 
For the remaining 20 PNe, the average iV c derived from various 
optical CEL diagnostics were used. For individual ionic species, 
we adopt the same effective reco mbination coeffici ents and assume 
the same case [Case A and B] as ILiu et al] d200fj|) . except for low 
I transitions of N II, for which the more recent calculations of 
iKisielius & Stor ey (2002) are used. For ionic species for which a 
host of lines have been detected, such as N II, O II and Ne II, the 
results derived from individ ual transit i ons w ere averaged, accord - 
ing to procedures outlined in lLiu et ID d 19951) and lLiu et alj d200d) . 
The results are tabulated in Tablesl8ll9l[T0land[TTl 

4.3 Elemental abundances 

Total elemental abundances were derived from ionic abundances 
(both from CELs and ORLs) using the sam e procedures and ion - 
ization correction factors (ICFs) adopted in lWesson et al.l d2005b . 
They are tabulated in Table|13l Ionic and elemental abundances as 
well as ICFs for individual elements, /(X), of all PNe are presented 
in Tables[8]|9l[T0l[n]and[l2] 

No [Fe IV] lines have been detected in any of the nebulae and 
no ICF formula is available to obtain Fe/H elemental abundance 
from the Fe 2+ /H+ ionic abundance alone. As such, we have ex- 
cluded Fe in our analysis. T he Mg II A4481 ORL has been detected 
in seven nebulae. Following lBarlowetal]d2003h . we assumed that 
the Mg II A4481 has an effective recombination coefficient equal 
to that of the C II A4267 line in calculating Mg 2+ /H+ from the 
observed strength of the A4481 line. In addition, considering the 
unusually wide ionization potential interval occupied by Mg 2+ , we 
assumed that Mg/H = Mg + /H + , i.e. no need of applying any ICF 
corrections. As a consequence, no Mg 2+ /H+ ionic abundances are 
tabulated separately. Instead, only the total elemental abundances 
of magnesium are given in Table [14] (see below). 

Total elemental abundances deduced for each nebula are pre- 
sented in Table [T4l For ease of comparison, we list abundances de- 
rived from both CELs and ORLs in the same Table. Also tabulated 
in the Table are average abundances of individual elements derived 
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Table 8. Ionic and elemental abundances of helium from ORLs. Ionic and elemental abundances of oxygen from CELs and from ORLs. 



He CELs ORLs 

Nebula He+/H+ He 2 +/H+ He/H O+/H+ 2 +/H+ /(O) O/H 2 +/H+ /(O) O/H 

average A4686 AA3726,29 AA4959, 5007 average 



He 2-118 


8.00e- 


-02 


8.00e- 


-03 


8.80e- 


-02 


7.91e- 


-06 


1.99e- 


-04 


1.065 


2.20e- 


-04 


3.44e- 


-04 


1.108 


3.81e- 


-04 


M2-4 


1.16e- 


-01 






1.16e- 


-01 


5.63e- 


-05 


4.76e- 


-04 


1.000 


5.33e- 


-04 


8.99e- 


-04 


1.118 


l.OOe- 


-03 


M2-6 


l.Ole- 


-01 






l.Ole- 


-01 


6.28e- 


-05 


2.37e- 


-04 


1.000 


3.00e- 


-04 


5.55e- 


-04 


1.264 


7.02e- 


-04 


M3-7 


1.22e- 


-01 


1.50e- 


-03 


1.24e- 


-01 


2.54e- 


-04 


4.41e- 


-04 


1.008 


7.00e- 


-04 


1.94e- 


-03 


1.588 


3.08e- 


-03 


M 1-20 


9.50e- 


-02 


4.00e- 


-05 


9.50e- 


-02 


2.30e- 


-05 


3.57e- 


-04 


1.000 


3.80e- 


-04 


4.95e- 


-04 


1.064 


5.27e- 


-04 


NGC 6439 


1.12e- 


-01 


2.10e- 


-02 


1.33e- 


-01 


5.82e- 


-05 


3.96e- 


-04 


1.121 


5.09e- 


-04 


2.44e- 


-03 


1.286 


3.14e- 


-03 


H 1-35 


1.02e- 


-01 


2.00e- 


-04 


1.02e- 


-01 


5.21e- 


-05 


2.78e- 


-04 


1.001 


3.31e- 


-04 


5.71e- 


-04 


1.188 


6.79e- 


-04 


M 1-29 


1.12e- 


-01 


3.30e- 


-02 


1.45e- 


-01 


1.08e- 


-04 


3.83e- 


-04 


1.187 


5.84e- 


-04 


1.13e- 


-03 


1.523 


1.72e- 


-03 


Cn 2-1 


l.lle- 


-01 


8.00e- 


-03 


1.19e- 


-01 


6.93e- 


-06 


4.79e- 


-04 


1.047 


5.09e- 


-04 


1.40e- 


-03 


1.062 


1.49e- 


-03 


H 1-41 


8.10e- 


-02 


2.10e- 


-02 


1.02e- 


-01 


1.99e- 


-05 


3.48e- 


-04 


1.166 


4.29e- 


-04 


1.78e- 


-03 


1.232 


2.20e- 


-03 


H 1-42 


1.09e- 


-01 


7.00e- 


-04 


l.lOe- 


-01 


3.21e- 


-05 


4.19e- 


-04 


1.004 


4.53e- 


-04 


9.64e- 


-04 


1.081 


1.04e- 


-03 


M2-23 


1.12e- 


-01 






1.12e- 


-01 


5.67e- 


-06 


2.58e- 


-04 


1.000 


2.64e- 


-04 


3.73e- 


-04 


1.021 


3.81e- 


-04 


M3-21 


1.14e- 


-01 


6.40e- 


-03 


1.20e- 


-01 


1.37e- 


-05 


6.30e- 


-04 


1.037 


6.68e- 


-04 


1.66e- 


-03 


1.059 


1.76e- 


-03 


H 1-50 


9.50e- 


-02 


l.lOe- 


-02 


1.06e- 


-01 


1.86e- 


-05 


4.28e- 


-04 


1.075 


4.80e- 


-04 


1.22e- 


-03 


1.122 


1.37e- 


-03 


M2-27 


1.27e- 


-01 


8.00e- 


-04 


1.28e- 


-01 


5.33e- 


-05 


6.86e- 


-04 


1.004 


7.42e- 


-04 


1.50e- 


-03 


1.082 


1.62e- 


-03 


H 1-54 


8.70e- 


-02 


3.00e- 


-05 


8.70e- 


-02 


5.43e- 


-05 


1.87e- 


-04 


1.000 


2.42e- 


-04 


4.77e- 


-04 


1.290 


6.15e- 


-04 


NGC 6565 


9.90e- 


-02 


1.50e- 


-02 


1.14e- 


-01 


1.50e- 


-04 


3.83e- 


-04 


1.098 


5.86e- 


-04 


6.47e- 


-04 


1.530 


9.90e- 


-04 


M2-31 


1.14e- 


-01 






1.14e- 


-01 


1.94e- 


-05 


4.37e- 


-04 


1.000 


4.57e- 


-04 












NGC 6567 


1.02e- 


-01 


9.00e- 


-04 


1.03e- 


-01 


1.24e- 


-05 


2.76e- 


-04 


1.005 


2.90e- 


-04 


6.04e- 


-04 


1.051 


6.35e- 


-04 


M2-33 


1.05e- 


-01 






1.05e- 


-01 


2.03e- 


-05 


4.93e- 


-04 


1.000 


5.13e- 


-04 


1.06e- 


-03 


1.041 


l.lOe- 


-03 


IC 4699 


8.00e- 


-02 


1.80e- 


-02 


9.80e- 


-02 


5.10e- 


-06 


2.67e- 


-04 


1.144 


3.11e- 


-04 


1.66e- 


-03 


1.166 


1.94e- 


-03 


M2-39 


1.12e- 


-01 






1.12e- 


-01 


6.55e- 


-05 


2.58e- 


-03 


1.000 


2.64e- 


-03 


9.14e- 


-04 


1.025 


9.37e- 


-04 


M2-42 


1.07e- 


-01 


3.00e- 


-04 


1.07e- 


-01 


3.21e- 


-05 


5.27e- 


-04 


1.001 


5.60e- 


-04 


l.lOe- 


-03 


1.062 


1.17e- 


-03 


NGC 6620 


l.lle- 


-01 


2.10e- 


-02 


1.32e- 


-01 


1.82e- 


-04 


5.11e- 


-04 


1.122 


7.78e- 


-04 


1.63e- 


-03 


1.521 


2.48e- 


-03 


VY 2-1 


1.29e- 


-01 


5.00e- 


-04 


1.30e- 


-01 


1.09e- 


-04 


5.50e- 


-04 


1.002 


6.60e- 


-04 


l.lle- 


-03 


1.200 


1.33e^ 


-03 


Cn 1-5 


1.25e- 


-01 


8.00e- 


-04 


1.26e- 


-01 


2.02e- 


-04 


4.89e- 


-04 


1.004 


6.93e- 


-04 


7.94e- 


-04 


1.418 


1.13e- 


-03 


M3-29 


l.OOe- 


-01 






l.OOe- 


-01 


5.88e- 


-05 


2.63e- 


-04 


1.000 


3.22e- 


-04 


5.76e- 


-04 


1.223 


7.05e- 


-04 


M3-32 


1.14e- 


-01 


1.05e- 


-02 


1.24e- 


-01 


2.58e- 


-05 


3.83e- 


-04 


1.060 


4.34e- 


-04 


6.80e- 


-03 


1.131 


7.70e- 


-03 


M 1-61 


1.04e- 


-01 






1.04e- 


-01 


3.66e- 


-05 


4.85e- 


-04 


1.000 


5.22e- 


-04 


9.49e- 


-04 


1.075 


1.02e- 


-03 


M3-33 


8.80e- 


-02 


1.70e- 


-02 


1.05e- 


-01 


6.87e- 


-06 


3.43e- 


-04 


1.125 


3.94e- 


-04 


2.25e~ 


-03 


1.147 


2.58e- 


-03 


IC 4846 


7.90e- 


-02 


2.00e- 


-04 


7.92e- 


-02 


5.01e- 


-06 


3.83e- 


-04 


1.001 


3.88e- 


-04 


5.91e- 


-04 


1.014 


6.00e- 


-04 



from the current sample of GBPNe, from other samples of GBPNe 
published in the literature, from the TLW (c.f. Section|7]( and KB94 
sample of GDPNe, as well as the solar photospheric abundances. 
Note that average abundances listed in the Table are logarithmic 
abundance averages (in units where H = 12), not averages of loga- 
rithmic abundances. 

4.4 Abundance uncertainties 

There are two major sources of uncertainty that affecting the de- 
rived element abundances: those propagating from line flux mea- 
surements and those arising from the analyzing methodology (e.g. 
those introduced via the application of empirical ICFs). Another 
potential source of error affecting in particular the current sample of 
objects is the uncertainty in the extinction law towards the Galactic 
center. The latter is difficult to characterize due to the lack of other 
independent means of measuring the extinction curve towards the 
Galactic center. 

Uncertainties in line fluxes will also affect abundance deter- 
minations via derivations of T c and N c . For most of our objects, 
we expect that flux measurements of strong lines (strengths larger 
than 0.008 relative to H/3) should be accurate to better than 5 
per cent. This encompasses essentially all diagnostic CELs, includ- 
ing the relatively weak [O III] A4363 and [N II] A5755 auroral 



lines, important for T c determinations. We estimated that uncer- 
tainties in plasma diagnostics of CELs should typically be no more 
than 4 per cent in T c and 30 per cent in 7Y e , and their total effect on 
deduced ionic abundances should be less than 6 per cent. There are 
a few exceptions where the abundances were based on low critical 
density far-IR fine-structure lines, which are more sensitive to A^ e 
variations. 

Emissivities of ORLs have only a weak dependence on T c and 
N B . However, those lines are generally much weaker than CELs 
and are often affected by serious line blending. For a few objects 
lacking high S/N, high resolution spectra, such as M 1-29, uncer- 
tainties in the derived ORL abundances of N, O and Ne could be as 
large as 30 per cent. The situation is much better for carbon, which 
has a few relatively strong well observed ORLs in the optical. Thus 
in most cases, we expect carbon abundances derived from ORLs 
to be accurate to better than 5 per cent. Comparable or even better 
accuracies are expected for helium abundances. 

Alexander & Balickl dl997h have investigated potential errors 
in abundance determinations that may result from applying the 
KB94's ICF scheme when only optical lines are observed, by con- 
structing a series of constant-density PN models. They found that 
the errors are smaller for small unresolved nebulae, and typically 
amount to < 10 per cent for helium and oxygen, < 25 per cent for 
nitrogen, sulfur and neon, but the abundances of argon are prob- 
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Table 9. Ionic and elemental abundances of carbon from CELs and from ORLs. 



CELs ORLs 

Nebula C 2 +/H+ C 3 +/H+ /(C) C/H C 2 +/H+ C 3 +/H+ /(C) C/H 

A1908 A1550 A4267 AA4187,4650 
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-05 


4.93e- 


-04 


1.289 


7.48e-04 


NGC 6565 


1.98e- 


-04 


1.530 


3.03e- 


-04 


3.27e- 


-04 


5.47e- 


-05 


1.392 


5.32e-04 


M2-31 












2.96e- 


-04 


3.38e- 


-03 


1.044 


3.84e-03 


NGC 6567 


7.80e- 


-04 


1.051 


8.20e- 


-04 


1.61e- 


-03 


7.00e- 


-03 


1.044 


9.00e-03 


M2-33 


3.94e- 


-04 


1.041 


4.10e- 


-04 


1.91e- 


-04 


l.lOe- 


-04 


1.041 


3.13e-04 


IC 4699 


1.02e- 


-04 


1.166 


1.19e- 


-04 


5.18e- 


-04 


6.82e- 


-05 


1.019 


5.97e-04 


M2-39 


6.09e- 


-03 


1.025 


6.25e- 


-03 


4.07e- 


-04 


5.91e- 


-04 


1.025 


1.02e-03 


M2-42 












1.58e- 


-04 


5.94e- 


-04 


1.060 


7.98e-04 


NGC 6620 


1.45e- 


-04 


1.521 


2.21e- 


-04 


8.25e- 


-04 


8.44e- 


-05 


1.355 


1.23e-03 


VY 2-1 


3.96e- 


-04 


1.200 


4.75e- 


-04 


3.06e- 


-04 


1.08e- 


-04 


1.197 


4.96e-04 


Cn 1-5 


1.05e- 


-03 


1.70e-04 1.412 


1.72e- 


-03 


1.49e- 


-03 






1.418 


2.10e-03 


M3-29 












3.33e- 


-04 






1.223 


4.08e-04 


M3-32 


2.61e- 


-04 


1.131 


2.95e- 


-04 


2.93e- 


-03 


2.30e- 


-04 


1.067 


3.37e-03 


M 1-61 












3.91e- 


-04 


1.79e- 


-04 


1.075 


6.13e-04 


M3-33 


9.96e- 


-05 


1.147 


1.14e- 


-04 


3.72e- 


-04 


2.62e- 


-04 


1.020 


6.47e-04 


IC 4846 


1.43e- 


-04 


1.014 


1.45e- 


-04 


l.lle- 


-04 


1.52e- 


-04 


1.013 


2.66e-04 



ably systematically overestimated. Given the large distance to the 
Galactic bulge, most PNe in our sample are quite compact, less than 
10" in diameter. A few objects in our sample have UV and/or IR 
spectra available. For those objects, we expect a better accuracy for 
nitrogen abundances. 

In summary, we estimate that the total error budgets of our cur- 
rent abundance determinations are: 0.12 dex for He/H from ORLs 
and O/H from CELs, 0.26 to 0.30 dex for N/H, S/H, Ne/H, Ar/H and 
Cl/H from CELs, 0. 12 dex for C/H from ORLs, 0.25 to 0.29 dex for 
O/H, N/H and Ne/H from ORLs. 



5 COMPARISONS WITH PREVIOUS STUDIES 

Hitherto several spectroscopic studies presenting approximately 
350 abundance determinations for more than 250 GBPNe have 
been published. In this section, we will compare our re- 
sults with those from some of the most recent ones, in- 
19881). 
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RPDM97 (Ratagetal. 
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SRM98 (Stasihs ka et al] 19981), C MKAS00 ( jCursinier et alj 20QI 
EC01 jEscudero & Costj|200ll) . ECM04 tescudero et al]|2004l) 
and EBW04 JExter et alJI 20041) . We will compare the general pat- 
terns of abundance obtained in Section IBTTI as well as individual 
measurements for common objects in Section |5T2l Unless other- 
wise specified, the current sample of 25 GBPNe, WL07, refers to 



the 23 GBPNe presented in the current paper (c.f. the first part of 
Table 1 141 with two peculiar objects, M2-23 and M 2-39 exclude d) 
plus M 1-42 and M 2-36 previously analyzed by Liu et al.1 (2001). 



Except for iLiu et al. (2001) presenting heavy element abun- 
dances deduced both from CELs and from ORLs for two GBPNe, 
essentially all abundance determinations of heavy elements for 
GBPNe published in the literature have been based on strong CELs 
only. Thus, comparisons will only be made for abundances derived 
from CELs. Also, carbon abundances have previously been mea- 
sured only for a very small number of GBPNe, so the comparisons 
will concentrate on helium, oxygen, nitrogen and neon. 



Note that much of the earlier work on GBPNe used the stan- 
dard extinction law of the general diffuse ISM, which is known 
to be inapplicabl e for the Galactic centre (e.g. UdalskJ 120031 : 
lRuffleetal]|2004l) . in particular in the UV. However, for analy- 
ses based on optical or infrared observations only, the uncertainties 
thus introduced were likely to be small. The recombination contri- 
bution to the [N II] A5755 auroral line was also uncorrected, lead- 
ing to overestimated T C ([N II]) and consequently underestimated 
N + /H + , + /H + and S + /H + ionic abundance ratios in some nebu- 
lae. 
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Table 10. Ionic and elemental abundances of nitrogen from CELs and from ORLs. 



CELs ORLs 



Nebula 


N+/H+ 
AA6548, 6584 


N 2 +/H+ 
A1750 


/(N) 


N/H 


N 2 +/H+ 
average 


N 3 +/H+ 
A4379 


/(N) 


N/I 


He 2-1 18 


1.54e-06 




27.85 


4.29e-05 


1.56e-04 




1.108 


1.73e 


M 2-4 


1.82e-05 




9.464 


1.72e-04 


3.84e-04 




1.118 


4.29e 


M 2-6 


1.12e-05 




4.782 


5.36e-05 










M3-7 


3.29e-05 




2.759 


9.08e-05 


1.07e-03 




1.588 


1.70e 


M 1-20 


5.29e-06 




16.51 


8.74e-05 


4.31e-04 


8.58e-05 


1.064 


5.50e 


NGC 6439 


3.14e-05 


3.38e-04 


0.304 


4.72e-04 


6.92e-04 


2.11e-04 


0.092 


9.67e 


H 1-35 


8.16e-06 




6.346 


5.18e-05 


2.46e-04 


1.06e-04 


1.187 


4.18e 


M 1-29 


7.34e-05 




5.393 


3.96e-04 


6.82e-04 


2.92e-04 


1.227 


1.20e- 


Cn2-1 


4,28e-06 




73.43 


3.14e-04 


2.98e-04 


4.92e-04 


1.013 


8.01e 


H 1-41 


3.50e-06 




21.53 


7.55e-05 


6.62e-04 


1.06e-04 


1.048 


8.05e 


H 1-42 


4.94e-06 




14.12 


6.98e-05 


6.88e-04 




1.081 


7.44e 


M2-23 


2.18e-06 




46.57 


1.02e-04 


7.81e-04 


6.05e-05 


1.021 


8.60e 


M 3-21 


8.25e-06 


8.38e-05 


0.420 


1.27e-04 


2.14e-04 


8.99e-05 


0.098 


3.25e 


H 1-50 


7.62e-06 




25.84 


1.97e-04 




8.92e-05 






M2-27 


4.56e-05 




13.91 


6.34e-04 


1.50e-03 




1.082 


1.62e 


H 1-54 


8.39e-06 




4.450 


3.73e-05 


4.15e-04 


1.22e-05 


1.289 


5.51e 


NGC 6565 


6.98e-05 




3.895 


2.72e-04 


9.03e-04 


3.87e-05 


1.345 


1.27e 


M2-31 


9.61e-06 




23.49 


2.26e-04 










NGC 6567 


2.20e-06 


3.51e-04 


0.024 


3.62e-04 


3.59e-04 


8.94e-06 


0.006 


3.70e 


M2-33 


2.23e-06 




25.25 


5.63e-05 


4.96e-04 




1.041 


5.16e 


IC 4699 


3.57e-07 




60.98 


2.18e-05 


2.32e-04 


1.05e-04 


1.016 


3.43e 


M2-39 


1.18e-05 




40.28 


4.77e-04 


7.84e-04 




1.025 


8.04e 


M2-42 


1.03e-05 




17.45 


1.80e-04 


4.39e-04 




1.062 


4.67e 


NGC 6620 


8.39e-05 




4.278 


3.59e-04 


6.41e-04 


1.47e-04 


1.305 


1.03e 


Vy 2-1 


2.68e-05 


2.08e-04 


1.002 


2.36e-04 


5.51e-04 




1.200 


6.61e 


Cn 1-5 


1.32e-04 


7.29e-04 a 


1.004 


8.65e-04 


4.93e-03 




1.418 


6.99e 


M3-29 


1.74e-05 




5.468 


9.50e-05 










M 3-32 


5.27e-06 




16.83 


8.87e-05 


2.24e-03 


3.04e-04 


1.063 


2.70e 


M 1-61 


8.97e-06 




14.23 


1.28e-04 


4.70e-04 


8.51e-05 


1.075 


5.97e 


M3-33 


4.19e-07 




57.27 


2.40e-05 




9.80e-05 






IC 4846 


1.58e-06 




77.45 


1.22e-04 


1.25e-04 




1.014 


1.27e 



a Derived from the [N III] 57/jm infrared fine-structure line. 

5.1 Comparisons of general patterns of abundances 

W88 presented elemental abundances for a sample of 49 GBPNe, 
amongst which five are in common with our sample: M 2-27, M 2- 
33, M 3-7, M 3-33 and Cn 2-1. Her sample yields an average oxy- 
gen abundance of 8.73, very similar to the value of 8.70 deduced 
from the current sample (c.f. Table I14t , For helium, her sample 
yields an average abundance of 11.06 (11.05 for non-Type I PNe), 
identical to the value deduced from the current sample. Note how- 
ever, collisional excitation of He I was not corrected for in W88 - 
the effects typically amount to a few per cent. She did not present 
nitrogen abundances separately. From the N/O abundance ratios 
listed in her paper, we derive an average N/O ratio of 0.24, much 
less than the average value of 0.40 deduced from the current sam- 
ple, suggesting that nitrogen abundances obtained by W88 are sys- 
tematically lower than ours by almost a factor of two. 

RPDM97 presented abundances for a sample of about 110 
GBPNe, deduced from their own observations or from line fluxes 
published in the literature. They used a grid of photoionization 
models to estimate ICFs, which were then used to calculate ele- 
mental abundances from abundances of observed ionic species. The 
average abundances calculated from values of individual nebulae 
presented in Table 5 of RPDM97 are listed in Table [H For He, 
O and S, their results are in excellent agreement with ours (within 
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0.05 dex). For N, CI and Ar, the differences are somewhat larger. 
For CI and Ar, measurement uncertainties may play a role - weak 
lines of [CI III] and [Ar IV] might not be well measured in their 
low-resolution spectra. For nitrogen, parts of the discrepancy may 
arise from uncertainties in ICFs. For CEL abundance analyses, only 
the trace ionic species N + is seen in the optical wavelength range. 
More abundant N ++ is only observable in the UV or the far in- 
frared. Thus ICF corrections for nitrogen are often quite large and 
uncertain when only optical data are available. 

SRM98 c ombin ed spectroscopic data published by 
lAller & Kevesl Jl987l) . W88 and RPDM97 and selected 90 
objects with a 6 cm radio flux density less than 100 mjy, con- 
sidered to belong to the Bulge population at the 90-95 per cent 
confidence level. They used an empirical abundance determination 
scheme similar to our CEL analyses and adopted the same set 
of ICFs given by KB94. However, they applied T c derived from 
either the [O IIl] n a or the [N Il] na ratios for all ionic species in 
abundance calculations. For their GBPNe sample containing 85 
objects with oxygen abundance measured, they derived a mean 
oxygen abundance 8.48 with a standard deviation (a) of 0.434. For 
32 brighter objects in their sample, the average abundance is 8.67 
(<t = 0.213). The latter is in good agreement with our result (c.f. 
Table [Hi. 



Elemental abundances of GBPNe from ORLs 

Table 11. Ionic and elemental abundances of neon from CELs and from ORLs. Ionic and elemental abundances of argon from CELs. 
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CELs ORLs CELs 

Nebula Ne 2 +/H+ /(Ne) Ne/H Ne 2 +/H+ Ne/H Ar 2 +/H+ Ar 3 +/H+ Ar 4 +/H+ /(Ar) Ar/H 

A3868 average A7136 AA471 1,4740 A7006 



He 2-118 


4.46e- 


-05 


1.108 


4.94e- 


-05 


1.55e- 


-03 


1.72e- 


-03 


4.29e- 


-07 


7.64e- 


-08 






1 


.037 


5.24e- 


-07 


M 2-4 


1.22e- 


-04 


1.118 


1.37e- 


-04 










2.56e- 


-06 


1.07e- 


-07 






1 


.118 


2.98e- 


-06 


M 2-6 


4.30e- 


-05 


1.264 


5.44e- 


-05 










6.42e- 


-07 


3.10e- 


-08 






1 


.264 


8.50e- 


-07 


M 3-7 


6.77e- 


-05 


1.588 


1.08e- 


-04 










6.33e- 


-07 










1 


.568 


9.93e- 


-07 


M 1-20 


6.04e- 


-05 


1.064 


6.43e- 


-05 










5.87e- 


-07 


7.33e- 


-08 






1 


.064 


7.03e- 


-07 


NGC 6439 


1.26e- 


-04 


1.286 


1.61e- 


-04 


1.35e- 


-03 


1.74e- 


-03 


2.27e- 


-06 


8.06e- 


-07 


l.Ole- 


-07 


1 


.071 


3.40e- 


-06 


H 1-35 


3.50e- 


-05 


1.188 


4.16e- 


-05 


8.57e- 


-05 


1.02e- 


-04 


1.35e- 


-06 


1.88e- 


-08 






1 


.187 


1.63e- 


-06 


M 1-29 


9.39e- 


-05 


1.523 


1.43e- 


-04 










1.95e- 


-06 


8.96e- 


-07 






1 


.227 


3.49e- 


-06 


Cn 2-1 


9.46e- 


-05 


1.062 


l.Ole- 


-04 


4.99e- 


-04 


5.30e- 


-04 


1.24e- 


-06 


4.22e- 


-07 






1 


.013 


1.68e- 


-06 


H 1-41 


8.30e- 


-05 


1.232 


1.02e- 


-04 










7.26e- 


-07 


4.27e- 


-07 






1 


.048 


1.21e- 


-06 


H 1-42 


8.56e- 


-05 


1.081 


9.25e- 


-05 


1.24e- 


-03 


1.34e- 


-03 


9.03e- 


-07 


2.61e- 


-07 






1 


.076 


1.25e- 


-06 


M 2-23 


3.45e- 


-05 


1.021 


3.52e- 


-05 


1.45e- 


-04 


1.48e- 


-04 


6.44e- 


-07 


7.62e- 


-08 






1 


.021 


7.36e- 


-07 


M 3-21 


1.72e- 


-04 


1.059 


1.83e- 


-04 


8.65e- 


-04 


9.17e- 


-04 


1.81e- 


-06 


8.56e- 


-07 


4.80e- 


-08 


1 


.069 


2.90e- 


-06 


H 1-50 


l.lOe- 


-04 


1.122 


1.23e- 


-04 










9.87e- 


-07 


7.10e- 


-07 


5.69e- 


-08 


1 


.040 


1.82e- 


-06 


M2-27 


1.90e- 


-04 


1.082 


2.06e- 


-04 










3.08e- 


-06 


4.61e- 


-07 






1 


.077 


3.81e- 


-06 


H 1-54 


2.59e- 


-05 


1.290 


3.34e- 


-05 


1.86e- 


-04 


2.40e- 


-04 


6.45e- 


-07 


1.17e- 


-08 






1 


.289 


8.47e- 


-07 


NGC 6565 


1.09e- 


-04 


1.530 


1.67e- 


-04 


8.20e- 


-04 


1.26e- 


-03 


1.74e- 


-06 


4.04e- 


-07 


7.03e- 


-08 


1 


.345 


2.99e- 


-06 


M2-31 


9.37e- 


-05 


1.044 


9.79e- 


-05 










1.42e- 


-06 


3.68e- 


-07 






1 


.044 


1.86e- 


-06 


NGC 6567 


4.71e- 


-05 


1.051 


4.94e- 


-05 


2.56e- 


-04 


2.69e- 


-04 


3.74e- 


-07 


1.20e- 


-07 






1 


.006 


4.96e- 


-07 


M2-33 


9.88e- 


-05 


1.041 


1.03e- 


-04 


7.04e- 


-04 


7.33e- 


-04 


1.20e- 


-06 


6.91e- 


-08 






1 


.041 


1.32e- 


-06 


IC 4699 


5.32e- 


-05 


1.166 


6.20e- 


-05 


5.48e- 


-04 


6.39e- 


-04 


1.46b- 


-06 


4.10e- 


-07 






1 


.016 


1.90e- 


-06 


M2-39 


2.10e- 


-04 


1.025 


2.15e- 


-04 










3.01e- 


-06 


1.46e- 


-06 






1 


.025 


4.58e- 


-06 


M2-42 


1.19e- 


-04 


1.062 


1.27e- 


-04 










1.44e- 


-06 


1.59e- 


-07 






1 


.060 


1.70e- 


-06 


NGC 6620 


1.44e- 


-04 


1.521 


2.19e- 


-04 


1.26e- 


-03 


1.92e- 


-03 


2.97e- 


-06 


6.23e- 


-07 


1.16e- 


-07 


1 


.305 


4.84e- 


-06 


VY 2-1 


1.41e- 


-04 


1.200 


1.70e- 


-04 


4.04e- 


-04 


4.85e- 


-04 


2.41e- 


-06 


1.08e- 


-07 






1 


.128 


2.84e- 


-06 


Cn 1-5 


1.59e- 


-04 


1.418 


2.26e- 


-04 


5.06e- 


-04 


7.18e- 


-04 


2.34e- 


-06 


7.19e- 


-08 






1 


.181 


2.85e- 


-06 


M3-29 


5.81e- 


-05 


1.223 


7.10e- 


-05 










6.33e- 


-07 










1 


.223 


7.74e- 


-07 


M3-32 


1.13e- 


-04 


1.131 


1.28e- 


-04 


2.00e- 


-03 


2.26e- 


-03 


1.31e- 


-06 


4.12e- 


-07 






1 


.063 


1.83e- 


-06 


M 1-61 


1.21e- 


-04 


1.075 


1.31e- 


-04 


1.49e- 


-03 


1.60e- 


-03 


1.43e- 


-06 


1.06e- 


-07 






1 


.075 


1.65e- 


-06 


M3-33 


8.87e- 


-05 


1.147 


1.02e- 


-04 










6.17e- 


-07 


5.49e- 


-07 






1 


.017 


1.19e- 


-06 


IC 4846 


5.74e- 


-05 


1.014 


5.83e- 


-05 










8.0()e- 


-07 


2.29e- 


-07 






1 


.013 


1.04e- 


-06 



CMKAS00 presented spectrophotometric observations and 
determined element abundances for a sample of 30 GBPNe. T c 's 
derived from the [N Il] n a and [O IIl] lla ratios were adopted for 
low- and high-ionization zones, respectively, and TVe's were de- 
rived f rom the [S Il] lln ratio. They used ICFs given by Kopp en et al.1 
Jl99lh . which are the same as ours for He, N and Ne. Mean abun- 
dances of He, O, N, S, CI and Ar, calculated from their Table 3, 
are presented in Table [74l For N, O and Ar, their results are al- 
most identical to ours, with differences less than 0.05 dex. For S 
and CI, the differences are larger, but still within the uncertainties 
(~ 0.2 dex). For helium their sample yields an average abundance 
of 10.98, 0.08 dex lower than our value, and is the lowest amongst 
all GBPNe s amples listed in Tab le[74l CMKAS00 used the criteria 
identified bv lAckeretai] ( ll99lh to select GBPNe, being similar to 
ours. There are however only two objects in common between the 
two samples. 

EC01 and ECM04 present element abundances for samples 
of 45 and 57 GBPNe, respectively. Unlike other studies, EC01 and 
ECM04 determined ionic abundances from measured line s t rength s 
using analytic functions deduced by I Alexander & BalickJ 119971) . 
rather than by solving level populations of th e emitting ions. The 
only e xception was Cl ++ , not dealt with by lAlexander & Balick 
( 1997); for this they solved statistical balance equations to obtain 
emissivities of the [CI III] AA55 17,5537 lines as a function of T c 



and iV c . Average abundances from these two samples of GBPNe 
were calculated from Table 6 of EC01 and Table 4 of ECM04 and 
tabulated in Table[74l 

The sample of ECO 1 seems to yield much lower average abun- 
dances for nitrogen and oxygen, by respectively as much as 0.60 
and 0.23 dex compared to our sample and that of ECM04, too large 
to accounted for by observational uncertainties alone. The large 
discrepancies are difficult to understand given that nebulae in the 
EC01 sample were observed and analyzed in the same way as those 
in the ECM04 sample, and the latter yields average abundances 
which are in excellent agreement with ours for all elements, with 
the possible exception of Ar for which our value is 0.39 dex loweiQ. 
It is not clear to us whether there were any systematic errors in the 
analysis of EC01. There are no objects in common between the 
EC01 and other samples, including ours, and EC01 did not publish 
line fluxes. Hence it is difficult to check whether the differences are 
real, caused by e.g. a different population being sampled by EC01 
or are possibly artifacts. 

Finally in a recent study, EBW04 determined abundances for 



4 If we exclude the 4 PNe in the ECM04 sample that show an Ar abundance 
either lower than 5 or higher than 7.6, which seems unrealistic, then the 
difference in the average Ar abundances deduced from the two samples 
drops to 0. 1 1 dex 
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Table 12. Ionic and elemental abundances of sulphur and chlorine from CELs. 



CELs CELs 

Nebula S+/H+ S 2 +/H+ /(S) S/H Cl 2 +/H+ /(CI) Cl/H 

AA6716,6731 A6312 AA5517, 5537 



He 2-118 


l.lle- 


-07 


1.77e- 


-07 


2.127 


6.13e- 


-07 


2.71e- 


-08 


3.461 


9.37e- 


-08 


M 2-4 


5.87e- 


-07 


1.92e- 


-05 


1.520 


3.00e- 


-05 


1.69e- 


-07 


1.566 


2.65e- 


-07 


M 2-6 


2.76e- 


-07 


3.29e- 


-06 


1.255 


4.48e- 


-06 


6.35e- 


-08 


1.360 


8.64e- 


-08 


M 3-7 


7.24e- 


-07 


8.09e- 


-06 


1.105 


9.74e- 


-06 


1.64e- 


-07 


1.204 


1.98e- 


-07 


M 1-20 


1.68e- 


-07 


2.48e- 


-06 


1.802 


4.77e- 


-06 


4.46e- 


-08 


1.924 


8.58e- 


-08 


NGC 6439 


1.12e- 


-06 


6.14e- 


-06 


1.485 


1.08e- 


-05 


1.28e- 


-07 


1.755 


2.25e- 


-07 


H 1-35 


2.46e- 


-07 


6.94e- 


-06 


1.354 


9.74e- 


-06 


9.25e- 


-08 


1.402 


1.30e- 


-07 


M 1-29 


2.39e- 


-06 


6.64e- 


-06 


1.295 


1.17e- 


-05 


1.15e- 


-07 


1.762 


2.03e- 


-07 


Cn 2-1 


1.80e- 


-07 


4.18e- 


-06 


2.916 


1.27e- 


-05 


6.95e- 


-08 


3.042 


2.11e- 


-07 


H 1-41 


1.43e- 


-07 


2.67e- 


-06 


1.959 


5.52e- 


-06 


6.76e- 


-08 


2.064 


1.40e- 


-07 


H 1-42 


2.90e- 


-07 


3.53e- 


-06 


1.716 


6.57e- 


-06 


8.18e- 


-08 


1.857 


1.52e- 


-07 


M 2-23 


1.16e- 


-07 


2.40e- 


-06 


2.512 


6.33e- 


-06 


7.82e- 


-08 


2.633 


2.06e- 


-07 


M 3-21 


3.29e- 


-07 


6.00e- 


-06 


2.551 


1.62e- 


-05 


1.22e- 


-07 


2.691 


3.27e- 


-07 


H 1-50 


3.57e- 


-07 


3.51e- 


-06 


2.076 


8.02e- 


-06 


7.26e- 


-08 


2.288 


1.66e- 


-07 


M2-27 


1.04e- 


-06 


1.04e- 


-05 


1.708 


1.95e- 


-05 


2.01e- 


-07 


1.879 


3.78e- 


-07 


H 1-54 


1.93e- 


-07 


4.00e- 


-06 


1.232 


5.17e- 


-06 


6.62e- 


-08 


1.292 


8.55e- 


-08 


NGC 6565 


2.52e- 


-06 


7.01e- 


-06 


1.192 


1.14e- 


-05 


1.12e- 


-07 


1.622 


1.81e- 


-07 


M2-31 


5.02e- 


-07 


6.05e- 


-06 


2.014 


1.32e- 


-05 


l.Ole- 


-07 


2.181 


2.21e- 


-07 


NGC 6567 


7.89e- 


-08 


1.36e- 


-06 


2.012 


2.90e- 


-06 


3.48e- 


-08 


2.129 


7.40e- 


-08 


M2-33 


5.95e- 


-08 


4.70e- 


-06 


2.061 


9.80e- 


-06 


1.03e- 


-07 


2.087 


2.15e^ 


-07 


IC 4699 


3.02e- 


-08 


7.59e- 


-07 


2.744 


2.16e- 


-06 


3.96e- 


-08 


2.853 


1.13e- 


-07 


M2-39 


3.81e- 


-07 


4.69e- 


-06 


2.396 


1.21e- 


-05 


4.43e- 


-06 


2.591 


1.15e- 


-05 


M2-42 


4.96e- 


-07 


6.46e- 


-06 


1.833 


1.28e- 


-05 


1.35e- 


-07 


1.974 


2.70e- 


-07 


NGC 6620 


3.25e- 


-06 


9.40e- 


-06 


1.220 


1.54e- 


-05 


1.82e- 


-07 


1.641 


2.98e- 


-07 


VY 2-1 


8.88e- 


-07 


l.OOe- 


-05 


1.339 


1.46e- 


-05 


1.53e- 


-07 


1.457 


2.23e- 


-07 


Cn 1-5 


2.55e- 


-06 


8.35e- 


-06 


1.158 


1.26e- 


-05 


1.59e- 


-07 


1.511 


2.40e- 


-07 


M3-29 


3.30e- 


-07 


3.48e- 


-06 


1.300 


4.96e- 


-06 


8.38e- 


-08 


1.423 


1.19e- 


-07 


M3-32 


2.46e- 


-07 


3.23e- 


-06 


1.813 


6.29e- 


-06 


8.09e- 


-08 


1.951 


1.58e- 


-07 


M 1-61 


2.82e- 


-07 


2.65e- 


-06 


1.720 


5.04e- 


-06 


9.30e- 


-08 


1.904 


1.77e- 


-07 


M3-33 


1.93e- 


-08 


1.21e- 


-06 


2.688 


3.30e- 


-06 


4.30e- 


-08 


2.731 


1.17e- 


-07 


IC 4846 


1.09e- 


-07 


3.36e- 


-06 


2.968 


1.03e- 


-05 


7.09e- 


-08 


3.064 


2.17e- 


-07 



a sample of 45 GBPNe, using a reduction procedure, set of atomic 
parameters and an ICF scheme that are almost identical to ours. Av- 
erage abundances derived from their sample are in excellent agree- 
ment with ours, with differences less than 0.04 dex for He, O, S, 
about 0.1 dex for N and Ne and about 0.26 dex for Ar (c.f. Ta- 
ble dl. 



5.2 Comparisons of abundances for common objects 

Fig.[8]compares abundances of He, N and O determined in the cur- 
rent work with those published in the literature for common ob- 
jects that were also observed and analyzed in the five other sam- 
ples of GBPNe discussed above, but excluding the sample of ECO 1. 
CMKAS00 have only two objects in common with us, and the He 
and N abundances they obtained are about 0.03 and 0.21 dex lower 
than ours, whereas those of oxygen are 0.10 dex higher. The level 
of agreement is quite good, considering the measurement uncer- 
tainties involved. The numbers of common objects are 6, 20, 5 and 
10 for the W88, RPDM97, ECM04 and EBW04 samples, respec- 
tively. 

In term of systematic errors, abundances obtained by EBW04 
are in best agreement with out results, although the scatter is also 
relatively large. The average abundance differences between the 
two samples are 0.00 ±0.07, -0.09 ±0.13 and 0.10 ±0.46 dex for 



He, O and N, respectively, where the errors are standard deviations. 
For ECM04, the abundances they deduced for He and O for the five 
objects in common are in excellent agreement with our results, with 
average differences of 0.00 ± 0.10 and —0.02 ± 0.07 dex, respec- 
tively. However, our nitrogen abundances are significantly lower 
than their results, by —0.37 ± 0.21 dex on average. 

For objects in common with the W88 and with the RPDM97 
samples, abundances of oxygen and nitrogen obtained in those 
studies are in good agreement with our results. For the W88 sam- 
ple, the average differences are 0.03 ± 0. 10 and 0.26 ± 0. 19 dex for 
oxygen and nitrogen, respectively. For the RPDM97 samples, the 
corresponding values are 0.01±0.22 and 0.14±0.23 dex. However, 
both W88 and RPDM97 yield helium abundances that are system- 
atically lower than our values, by 0.05 ± 0.03 and 0.04 ± 0.04 dex 
on average, respectively. The differences are found to mainly arise 
from measured line fluxes - for the same object, helium line fluxes 
measured by W88 and by RPDM97 are systematically lower than 
our values. For objects in common with RPDM97, if we recal- 
culate helium abundances with our own method but using their 
line fluxes, the results are almost identical to their values. For the 
W88 sample, if we apply our analysis to her line data, the resul- 
tant helium abundances decrease by approximately 0.03 dex. This 
is caused by the fact W88 did not make any correction for the ef- 
fects of collisional excitation of He I lines. Both W88 and RPDM97 
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Table 13. Ionization correction factors (ICFs). 
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(cjtf)cels 
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73 
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applied non-linearity corrections to their line fluxes, measured, re- 
spectively, with the 3.9m Anglo-Australian Telescope using an Im- 
age Dissector Scanner as detector and with the ESO 3.6m or 1.52m 
telescopes also using an IDS. It is quite possible that both W88 and 
RPDM97 have over-corrected the effects of non-linearity, leading 
to underestimated fluxes for lines of low or intermediate strengths. 
If this was indeed the case, then intensities of the [O III] A4363 and 
[N II] A5754 auroral lines were also likely to have been underesti- 
mated, leading to underestimated T c , with apparent consequences 
to ionic abundances deduced for heavy elements. 



6 CELS VERSUS ORLS 

As described in SectionfJJ plasma diagnostic analyses and elemen- 
tal abundance determinations of PNe and H II regions using CELs 
and ORLs/continua have ubiquitously showed dichotomy between 
results deduced from those two types of emission excited by dif- 
ferent mechanisms. The discrepancies are generally characterized 
by AT C , defined as T c ([0 III])-T C (BJ), the difference between T e 
deduced from the [O III] collisionally excited lines and from the 
B aimer discontinuity of H I recombination spectrum, and by abun- 
dance discrepancy factor, adf, defined as the ratio of ionic abun- 
dances X 1+ /H + derived from ORLs and from CELs. Deep op- 
tical spectroscopic surveys for dozens of GDPNe and a number 
of Galactic and extragalactic H II regions have shown that values 
of T e (BJ) are systematically lower than those of T c ([0 III]) de- 
duced for the same objects and that adf's are always larger than 
unity, falling typically in the range 1.6-3, but with a significant 
tail extending to higher values, reaching ~ 70 in the most e xtrem e 
case. Except for the two GBPNe analyzed by iLiu et al. 



ORL studies of PNe have so far concentrated on relatively bright 
GDPNe. The sample of GBPNe presented in the current work thus 
provides an interesting environment to probe the behaviours of the 
aforementioned fundamental discrepancies for a different stellar 
population. In particular, previous studies of GDPNe show that rel- 
ative abundances of heavy elements, such as the C/O ratios deduced 
from ORLs are compatible to those derived from CELs. It is impor- 
tant to test whether the same result also holds for GBPNe. 

6.1 Electron temperatures 

In Figs. [10] and [TT] T C (BJ) is plotted against T c ([0 III]), AT C , 
T e (He I A7281/A6678) and against T e (0 II ORLs). Nebulae of 
sub-sample B, for which only low resolution spectra are available 
near the Balmer discontinuity region, do not behave obviously dif- 
ferently compared to objects of sub-sample A, although measure- 
ments for sub-sample B clearly suffer from larger uncertainties 
compared to the latter sub-sample. As a consequence, in analyses 
below we shall concentrate on patterns and relations revealed by 
sub-sample A. 

A majority of our sample GBPNe (15 out of 25) have T C (BJ) 
higher than T c [0 III] (c.f. Fig.[9j. For the sub-sample of high 
data quality, exactly half of the objects have T C (BJ) higher than 
TLdO III]) while the other half have lower T e (BJ). The result ap- 
pears to contradict to what previously found for GDPNe where for 
majority of them, T C (BJ) is systematically lower than T c ([0 III]). 
Note that the magnitude of Balmer jump is inversely proportional 
to T c , thus as T c increases, the Balmer jump becomes increasing 
difficult to measure. And because of the nonlinear dependence of 
T C (BJ) on the magnitude of Balmer jump, the resulting error of 
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Table 14. Total elemental abundances on a logarithmic scale where H = 12. 
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8.09 


9.03 


7.02 


5.29 


6.31 


7.72 
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11.09 
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a Excluded from statistical analyses (c.f. Section |3~3l . 

b Average abundances for 23 GBPNe of the current sample. 

c Average abundances for 23 GBPNe of the current sample plus M 1-42 and M 2-36 previously analyzed bv lLiu et"al] j200ll) . 

d Combined sample of 58 GDPNe (c.f. Section|7J. 

e lLodderdj2003l) . 
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Figure 9. T C (BJ) plotted against a) T o ([0 III]) and b) AT e = T e ([0 III])-T C (BJ) for GBPNe (filled circles and asterisks) and GDPNe (open circles). Filled 
circles and asterisks represent respectively GBPNe of sub-samples A and B, denned in Section |2TI In a), the dotted-dashed diagonal line represents identical 
T C (BJ) and T e ([0 III]). 



T C (BJ) for a given uncertainty of BJ is not symmetric - it biases 
towards higher values. It seems to us, at least parts, if not all, of 
the discrepant behaviour of GBPNe compared to GDPNe as seen 
in Fig.[9] are due to lower quality of the data available for GBPNe, 
which are generally fainter than nearby GDPNe. 

While the left panel of Fig.|9] shows no obvious correlation 
between T C (BJ) and T c ([0 III]), the right panel shows that there 
is a tight anti-correlation between T C (BJ) and AT a and different 
samples seem to yield similar results. Both panels however re- 
flect the fact that while T (BJ) varies over a wide range of val- 
ues, T e ([0 III]) falls wi thin a muc h narrow range, as previously 
observed in GDPNe (c.f. Liu 2006a b and references therein). 

Fig.fTOlcompares T C (BJ) with T e (He I A7281/A6678), T c de- 
rived from the ratio of He I singlet recombination lines A7281 and 
A6678. For all PNe plotted in the Figure, both disk and bulge ones, 
T e (He I A7281/A6678) is invariably lower than T C (BJ), in consis- 
te nt with the expec tations of the two-abundance model proposed 
bv lLiu et all d200d) . which predicts t hat T e (Q II ORLs) < T^He I 
OR Ls) < T e (H I BJ) < T c ([0 III]) dPequignot et al]|2002[ [20031 ; 
lLivj|2003r) . exactly opposite to the expectations of the paradigm of 
temper ature fluctuations w hich predict T c (He I A7281/A6678) > 
T C (BJ) dZhang et alj2005t) . 

Fig.Qjjplots T c (0 II ORLs) against T C (BJ), and again, in con- 
sistent with what previously found for GDPNe, values of T c (0 II 
ORLs) are considerably lower. 

About half PNe in our GBPNe sample have an [O II] 
A4089/A4649 line ratio higher than the low T c limit of 288 K, the 
lowest T e at which the effective re combination coefficients for O II 
are available. iTsamis et ai]( l2004h detected a faint emission line at 
41 16 A in three bright PNe, NGC 3242, NGC 5882 and NGC 6818, 
and attributed it to Si IV 4s 2 Si/ 2 -4p 2 Pi /2 A4116.1. They cor- 
rected for the contribution of the other component of this multiplet, 
the 4s 2 Si/2-4p 2 P 3/2 A4088.8 line, to the observed intensity of the 
O II A4089.3 line. After the corrections they obtained an T e (0 II 
ORLs) of 2600, 8700 and 2900 K for NGC 3242, NGC 5882 and 
NGC 6818, respectively. For comparison, the H I Balmerjump tem- 
peratures they obtained for the three PNe are 10,200, 7800 and 
12,100 K, respectively, and the corresponding values for T c ([0 III]) 
are 11,700, 9400 and 13,300 K. The 4116 A feature has not been 



detected in NGC 615 3 dLiu et alJl200Ct adf = 9.7), in bulge PNe 
M 1-42 and M2-36 dLiu et alj|200ll; adf = 22.0 and 6.9, respec- 
tively) and in H f 2-2 dLiu et alj|2006l ; adf = 71.2). For NGC 7009 
dLiu et al 1 19951 ; adf = 4.7), our unpublished deep spectra reveals a 
weak feature at 4116 A. And if we apply similar corrections, then 
T c (0 II ORLs) yielded by the A4089/A4649 ratio increases from 
420 K before the correc tion to 13 1 K. F or the sample of PNe ob- 
served and analyzed by iLiu et"aL ( 2004a. b), a line 4116.24 A rest 
wavelength was detected with a dereddened intensity of 0.029 (in 
units where H/3 = 100) in the spectrum of NGC 6884 and was 
identified as [Ni V] A41 17.20. If we instead attribute it to Si IV 
A41 16, then T c (0 II ORLs) yielded by the A4089/A4649 ratio will 
increase from the value of 3040 K before the correction to 22, 100 K 
after the corre ction. Finally, among st the large number of sample 
PNe studied by Wess on etall d2005h . the only PN that has a feature 
detected at 41 16 A is Vy 2-2, with an observed intensity of 0.178 in 
units of H/3 = 100 after reddening corrections. However, the fea- 
ture is unlikely to be Si IV A41 16 for two reasons: 1) If it is Si IV 
A41 16, then the other component of the multiplet at A4089 should 
have an intensity of 2 x 0.178 = 0.356, which is over a factor of 
three higher than the measured total intensity of the A4089 feature 
(0.098); 2) The ++ /H + ionic abundance ratio derived from the 
O II A4089 line, assuming no contamination from Si IV, is in good 
agreement with values deduced from other O II ORLs. The case 
of Vy2-2 suggests that the identification of the 41 16 A features 
detected in some PNe as the Si IV A4116 is not completely secure 
and further observations of better S/N's and spectral resolutions are 
needed to clarify the situation. 



6.2 Abundance discrepancy factors 

The abundance discrepancy factor (adf) for a heavy element 
ionic species X 1+ is defined as ratio of ionic abundance ra- 
tios X i+ /H+ derived from ORLs and from CELs, adf(X i+ ) = 
(X 1+ /H + )orls/(X 1+ /H + )cels- Previous studies of several sam- 
ples of GDPNe show that adf varies over a wide range of values, 
from near unity (i.e. the ionic abundance ratios derived from the 
two types of emission lines agree) to over a factor of ten, peaking 
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Figure 8. He, N and O abundances derived in the current work are compared 
to those in the literature for common objects. Filled circles: W88; Open cir- 
cles: RPDM97; Open squares: CMKAS00; Open triangles: ECM04; Aster- 
isks: EBW04. 
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as in Fig. [9] 



around a factor of two but with a tail extending to much higher val- 
ues, re aching a factor of ~ 70 in the most extreme case (c.f. iLivj 
l2006allbh . 

In Fig. [12] we plot and compare the accumulative distribu- 
tions of adf for disk PNe from the TLW (c.f. Section |7J sample 
and for bulge PNe from the current sample. The two distributions 
show some differences for adf(0 2+ ) between 0.5 and 0.8 dex. A 
Kolmogorov-Smirnov Two-sample test, with the Null Hypothesis 
that the two samples have the same distributions of adfs, yields 
a Kolmogorov-Smirnov statistic D of 0.1808 with a probability 
of the two distributions having a difference exceeding \D\ under 
the Null Hypothesis, P(> \D\) — 0.59, which is equal to 1 mi- 
nus the confidence level. In other words, the risk is high, nearly 
60 per cent, to reject the Null Hypothesis. The analysis thus shows 
that to a fairly good confidence level, both GDPNe and GBPNe 
seem to yield similar distributions of adf(0 2+ ). 

In Fig. [T3] we plot total elemental abundances of oxygen de- 
rived from CELs and from ORLs against adf(0 2+ ). The Figure 
shows that except for a few PNe with O/H less than 8.2, oxygen 
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Figure 12. Comparison of the accumulative distributions of log adf(0 2+ ) 
of GDPNe (solid line) and of GBPNe (dotted line). 
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Figure 13. Elemental abundances of oxygen derived from ORLs (triangles) 
and from CELs (circles) plotted against adf(0 2 + ) for GBPNe (filled sym- 
bols) and GDPNe (o pen symbols). The solar photospheric oxygen abun- 
dance jLoddersl2003l) is marked. 



abundances derived from CELs fall into a relatively narrow range, 
consistent with the general consensus that oxygen is neither cre- 
ated nor destroyed by any significant amount in LIMS. By con- 
trast, oxygen abundances yielded by ORLs, apart from being al- 
ways higher than the corresponding values derived from CELs for 
the same object, cover a wide range of values, reaching about a fac- 
tor of 25 higher than the solar value in the most extreme case. Such 
high abundances are extremely difficult to explain in terms of the 
current theory of evolution of low- and intermediate-mass (single) 
stars, and argue strongly against the paradigm of temperature fluc- 
tuations as the underlying physical cause of the large adf values, 
as the scenario implicitly assumes that the higher abundance values 
yielded by ORLs, being insensitive to t he presence of temperature 
fluctuations, are the correct values (c.f. Liu 2006b). Together with 
Fig. [9] Fig. [T3] also implies that whatever effects are causing the 
large values of adf observed in PNe, it seems they manifest them- 
selves mainly by affecting emission of ORLs/continua, leading to 
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Figure 14. adf(0 2 +) plotted against AT C . For the meaning of different 
symbols, c.f. Fig. [9] The solid line shows the linear fit obtained bv lLiu et"al] 
1200 ll) . 



large variations in T c and ionic abundances that are inferred from 
emission of ORLs/continua. 

Fig. 1141 plots log adf(0 2+ ) against AT . Also overplotted 
is the linear fit obtained for a much smaller sample of PNe by 
iLiu et all d200ll) who first noticed that the two quantities are pos- 
itively correlated. The result suggests that both phenomena, i.e. 
that T derived from the [O III] forbidden lines are systematically 
higher that those derived from the Balmer discontinuity of the H I 
recombination spectrum, and that ionic abundance ratios of heavy 
elements derived from ORLs are systematically higher than de- 
duced from CELs, are related and probably caused by the same 
underlying physical mechanism. Furthermore, from the discussion 
above, we may conclude that the mechanism, whatever it is, affects 
mainly the emission of ORL/continua, rather than that of CELs. 
The results fit nicel y with the expec tations of the two-abundance 
model proposed by ILiu et alj d2000h . but contradict those of the 
paradigm of temperature fluctuations which affect mainly the pro- 
duction of emission of CELs. 

The mean surface brightness reflects the evolutionary stage 
of a PN. As the PN evolv es and expands, its density drops and 
thus the surface brightness. Gar nett & Dinersteinl ( feOOll) found that 
adf(0 2+ ) tends to be larger for big ger, lower surfac e brightness 
PNe. A similar trend was found bv lLiuetalJ (2004b). In Fig. Q3] 
we plot adf(0 2+ ) against the nebular 6 cm radio surface brightness, 
defined as F(5GHz)/0 2 , where F(5GHz) is nebular 6 cm ra- 
dio flux density in mjy and 9 is angular diameter in arcsec (c.f. 
Table [T}. Using radio flux density instead of H/3 flux has the ad- 
vantage that the former is not affected by interstellar extinction. 
Although with large scatter, the diagram shows some evidence that 
adf increases with decreasing surface brightness, consistent with 
prev ious results. 

IZhang et alj (2004) found a positive correlation between AT e 
and helium abundance, in accordance with the expectations of the 
two-abundance model. Our current data for GBPNe exhibit a simi- 
lar trend, albeit with larger scatter due to the larger uncertainties 
associated with the current determinations of T C (BJ) (c.f. Q2A\ . 
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Figure 15. adf(0 2 +) plotted against radio band surface brightness at 6 cm. 
For the meaning of different symbols, c.f. Fig. [9] For a few PNe that have no 
measurements of 5 GHz radio flux density, we have adopted 1 .4 GHz fluxes 
and converted the values to 5 GHz fluxes assuming optically thin emission. 
Those objects are marked with triangles. 
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Figure 17. Comparison of the distributions of log adf for 2+ /H + (solid 
curve), C 2+ /H+ (dashed-dotted histogram), N 2 +/H+ (dotted histogram) 
and Ne 2+ /H+ (dashed histogram) ionic abundance ratios. 
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Figure 16. Helium abundance plotted against adf(0 2+ ). For the meaning of 
di fferent symbols, c .f. Fig. [9] The solid line represents the linear fit obtained 
bv lLiu et alj <2004bl) . 



Fig.QUlplots He abundance against adf(0 2+ ). T he solid line in the 
Figure shows the linear regression obtained by iLiu et al. (2004b) 
for disk PNe. Again, while the scatter is large, there is a discernible 
positive correlation between the two quantities. 

It has been shown that apart from oxygen, the large discrep- 
ancies between ionic abundances relative to hydrogen derived from 
ORLs and from CELs exist for the oth er three abundant second- 
row e lements, carbon, nitrogen and neon iLiu et alJl995l ; lLuo et al.l 
1200 lh . In addition, for a given PN, values of the adf deduced for 
the four most abun dant second-row heavy elements are compara- 
ble fciu et alj[2OO0h . In other words, both CEL and ORL analyses 



yield similar relative abundance ratios of heavy elements, such as 
C/O, N/O and Ne/O. Fig. [17] shows the distributions of log adf(X), 
where X = C 2+ /H + , N 2+ /H + , 2+ /H+ and Ne 2+ /H+, for all 
Galactic PNe for which a measurement of the quantity is avail- 
able. The distribution of adf(0 2+ ) peaks strongly near 0.3 dex, with 
a FWHM of 0.26 dex and an extended tail towards higher values. 
For C 2+ /H + and N 2+ /H + , the distributions peak around 0.5 dex 
with a FWHM of about 0.7 and 1.0, respectively. For Ne 2+ /H + , 
it peaks around 0.8 dex with a FWHM of 1.2 dex. In Fig. [18] 
we plot the adf deduced for 2+ /H + against those of C 2+ /H + , 
N 2+ /H + and Ne 2+ /H+. The Figure shows that adf(C 2+ ) follows 
closely adf(0 2+ ), whereas adf(N 2+ ) and adf(Ne 2+ ) are system- 
atically higher than adf(0 2+ ). Heavy element abundance ratios, 
C/O, N/O and Ne/O, derived from CELs and from ORLs, are com- 
pared in Fig. [19] In line with the patterns revealed in Fig.[T8] both 
CELs and ORLs yield compatible C/O abundance ratio, whereas 
for N/O and Ne/O, the ratios deduced from ORLs are systemat- 
ically higher, a direct consequence of higher values of adf found 
for N 2+ /H + and Ne 2+ /H + . The average C/O ratios found from 
ORLs are -0.11 ± 0.05 for GDPNe (39 objects) and -0.38 ± 0.06 
for GBPNe (13 objects). The corresponding results from CELs are 
-0.19 ± 0.05 and -0.35 ± 0.10 for GDPNe and GBPNe, re- 
spectively. Injhe scenario of two-abundance scenario proposed by 
ILiu et alj < |2000|) . ORLs and CELs trace respectively the cold, H- 
deficient inclusions and the hot, diffuse medium of normal com- 
position. The fact that both ORLs and CELs yield compatible C/O 
ratios suggests that the postulated H-deficient inclusions have the 
same C/O ratio as the diffuse medium and therefore unlikel y to 
be nucleo-processed material as proposed bv llben et al. Ul983l) for 
H-deficient knots observed in Abell 30 and Abell 78. Deep spec- 
troscopy and ORL abundan ce analysis of H-deficient knots in 
Abell 30 dWesson et al.ll2003l) also show that the knots are O-rich 
rather than C-rich, in contradiction to p redictions of the born-again 
scenario proposed bv llben et alj < ll983l) . 

Fig. [19] also shows the C/O ratio varies over a factor of ten, 
from —0.5 to +0.5 dex, and that C/O ratios of GBPNe are system- 
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Figure 19. Comparison of heavy element abundance ratios C/O, N/O and Ne/O derived from ORLs and from CELs, for GBPNe (filled circles) and GDPNe 
(open circles). The solid lines represent y = x. 



atically lower than those of GDPNe. This important result will be 
discussed further in the next Section. 

For the four abundant second-row elements, adf(0 2+ ) is prob- 
ably the best determined. First of all, high quality calculations of 
effective recombina tion coefficient s for the whole range of O II 
ORLs are available dLiu et alJI 19951) . Those calculations have been 
extended down to a T c a s low as 100 K and even to 288 K for 
some selective transitions dBastinl2006t) . Secondly, most prominent 
O II ORLs fall within the narrow 4000-5000 A optical wavelength 
range where one also finds the three [O III] CELs, the AA4959,5007 
nebular lines and the A4363 auroral line, upon which CEL analysis 
relies. All those lines can thus be measured simultaneously using 
a single instrument and sampling the same area of a given nebula. 
This avoids or minimizes errors caused by effects such as uncer- 
tainties in flux calibration and reddening c orrections. For c arbon, 
good atomic data also exist for C II ORLs jDavev et alfeood) and 
the lines are generally well measured, in particular the relatively 
strong 3d-4f A4267 transition, yielding fairly accurate C 2+ /H + ra- 
tios. On the other hand, for CEL analysis, the derivation of carbon 
ionic abundance ratios has to rely on UV CELs, most importantly, 
the C III] AA1907,1909 for C 2+ /H+ and C IV AA1548,1550 for 
C 3+ /H + . Thus except for a number of extended PNe for which the 
opti cal spectra have been obtained by s canning the whole nebula 
(e.g. iLiu et al. 2000; Tsamis et al . 2004), and for those nebulae of 
small angular sizes (< a few arcsec, such as most GBPNe) and 
nearly symmetric, the adf(C 2+ ) deduced may suffer from errors. 
These uncertainties result from the fact that the UV and optical ob- 
servations may not be sampling the same portions of the nebula as 
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Figure 20. Magnesium abundance, Mg/H m Mg 2 +/H + , with the later de- 
rived from the Mg II A4481 ORL, plotted against adf(0 2 +). The solar pho- 
tospheric magnesium abundance (Lodders 2003) is marked. For an expla- 
nation of different symbols, c.f. Fig.|9] 



well as from potential uncertainties in flux normalization and red- 
dening corrections. 

Good atomic data n ow exist fo r I ^ 2 transitions of N II 
dKisielius & Storevll2002l) and Ne II jKisielius et ai]| 19981) ORLs, 
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but not those arising from levels of higher orbital angular mo- 
menta, such as 3d~4f lines. N II and Ne II ORLs are generally much 
more difficult to observe than O II ORLs due to lower abundances 
of nitrogen and neon compared to oxygen. Most of the strongest 
Ne II 3-3 transitions fall at shorter wavelengths below 4000 A. 
Whereas Ne 2+ /H + can be deduced from strong optical [Ne III] 
AA3868,3967 nebular lines, N 2+ /H+ can only be obtained from 
either the relatively weak N III] A1751 intercombination line or 
from the [N III] 57^im far-infrared lines. All those factors make 
the current estimates of adf(N 2+ ) and adf(Ne 2+ ) subject to large 
uncertainties. We believe that a major fraction, if not all, of the 
larger values of adf found for N 2+ /H + and Ne 2+ /H + compared to 
adf(0 2+ ) (Fig.ll8t and consequently, higher N/O and Ne/O ratios 
yielded by ORLs compared to those obtained from CELs (Fig.|19ll, 
are caused by those uncertainties and are therefore artificial. Better 
observations and atomic data for N II and Ne II ORLs are clearly 
needed. More accurate determinations of N/O and Ne/O ratios us- 
ing ORLs are of great interest to unveil the nature and origins of 
the postulated H-deficient inclusions embedded in the diffuse neb- 
ula that are responsible for large values of adf observed in PNe. 

Magnesium is the only third-row element for which an ORL 
from one of its dominant ionic species has been detected, i.e. the 
Mg II 3d-4f A4481 line yieldi ng Mg 2+ /H+ To a good approx- 
imation, Mg/H ~ Mg 2+ /H+ ( iBarlow et all 120031) . Mg/H abun- 
dances thus deduced are now available for a total of 32 PNe, in- 
clu ding 24 GDPNe an d 8 GBPNe (excluding the peculiar M2- 
24; IZhang & Liull2003l) . The results are plotted in Fig. 1201 against 
adf(0 2+ ). Mg/H abundance is particularly difficult to obtain from 
CELs, as it has to rely on weak lines in the UV or in the infrared. 
Thus ad f(Mg 2+ ) cannot be deduced. However, as previously no- 
ticed bv lBarlow et al.l d2003h . in stark contrast to C, N, O and Ne, 
for a wide range of adf(0 2+ ), Mg/H deduced from the Mg II A448 1 
ORL fall in a narrow range compatible with the solar photospheric 
value. The 24 GDPNe plotted in Fig.[20]yield an average Mg/H ra- 
tio of 7.56 dex with a standard deviation of .22 dex, almost identi- 
cal to the solar value of 7.55 jLoddersll2003h . The 8 GBPNe yield 
a slightly higher average value, of 7.71 dex. We will discuss mag- 
nesium abundance in the context of the chemical evolution of the 
Galaxy in the next Section. 

6.3 Summary 

In summary, the current detailed plasma diagnostics and abundance 
analyses contrasting results obtained from CELs and those deduced 
from ORLs and continua for a large sample of GBPNe have yielded 
the same phenomena and strengthen the conclusions previously 
found for GDPNe: 

1) Nebular electron temperatures deduced from recombination 
lines or continua are systematically lower than those deduced from 
the classic [O III] nebular to auroral forbidden line ratio. In addi- 
tion, there is strong evidence that T c derived from He I ORL ratios, 
T c (He I), are systematically lower than T C (BJ), T c deduced from 
the H I recombination continuum Balmer discontinuity, and that T c 
given by recombination line ratios of heavy element ions, such as 
from O II ORL ratios, are even lower, yielding values that are just 
a few hundred K in some cases. 

2) Heavy element abundances relative to hydrogen derived 
from ORLs are always higher than values derived from CELs. This 
is true for all abundant second-row elements, C, N, O and Ne. For a 
given nebula, values of the adf found for those individual elements 
are comparable, leading to similar values of abundance ratios, C/O, 
N/O and Ne/O deduced from ORLs and from CELs. By contrast, 



no enhancement of ORL abundances relative to CEL values is ap- 
parent for magnesium, the only third-row element that has been 
studied using ORLs. There is also evidence that objects exhibiting 
particularly large values of adf also have higher helium abundances. 

3) The observed discrepancies between T c and heavy element 
abundances derived from ORLs and from CELs are correlated. Ob- 
jects showing large values of adf tend to have large T c discrepan- 
cies, suggesting that the two phenomena are related, and are likely 
caused by the same underlying physical process. 

4) The large discrepancies between T a derived from CELs 
and from ORLs appear to be mainly caused by abnormally low 
values yielded by recombination lines or continua. For the large 
number of nebulae analyzed, T a derived from CELs fall in a nar- 
row range around 10,000 K, consistent with what is predicted by 
photoionization models for nebulae of "normal" chemical composi- 
tion, i.e. that of Sun, a typical G-dwarf in the Galaxy. Similarly, the 
large discrepancies between heavy element abundances deduced 
from ORLs and from CELs are largely caused by the abnormally 
high values obtained from ORLs. Excluding a few nebulae that are 
known to belong to the halo population and thus have lower heavy 
element abundances, abundances of most PNe deduced from CELs 
fall in a narrow range of slightly sub-solar. By contrast, ORLs yield 
much higher heavy element abundances, reaching more than a fac- 
tor of ten higher than the solar values in some extreme cases. It 
appears that whatever is causing the large differences between re- 
sults obtained from ORLs and from CELs, the main effect of the 
process is to enhance emissivities of recombination lines and con- 
tinua, but hardly affect those of CELs. The results also suggest that 
heavy element abundances deduced from ORLs are unlikely to be 
real, or representative of the bulk composition of the whole nebula. 
This is in stark contradiction to the prediction of the paradigm of 
temperature fluctuations, which implicitly assumes that the higher 
abundances yielded by ORLs are the correct ones. 

5) All the above results point to the existence of a cool compo- 
nent of meta l-rich plasma em bedded in the diffuse nebula, as first 
proposed bv lLiu et all ( feood) . The cool inclusions are responsible 
for most of the observed fluxes of ORLs emitted by heavy element 
ions, yet contribute little to the fluxes of CELs. As the nebula ex- 
pands and ages, the diffuse emission fades and emission from the 
embedded cool inclusions becomes more apparent. 



7 DISCUSSION 

In this Section, we compare elemental abundances deduced for 
the current sample of GBPNe with those derived previously for 
GDPNe. We investigate the abundance patterns and discuss the re- 
sults in the context of the nucleosynthesis and dredge-up processes 
in LIMS and their roles in the chemical evolution of the Galaxy. 
Of the GBPNe observed in the current work (Table [T), we have 
excluded in our analysis the peculiar nebulae M 2-23 and M 2-39 
(c.f. §13. 3t. but i ncluded M2-36 and M 1-42 previously studied by 



iLiu et alj j200lh . The sample, WL07, consists of 25 bulge PNe in 
total. 

As described in §Q] deep spectroscopy, contrasting plasma 
diagnostic and abundance analysis results obtained from the tra- 
ditional method based on strong CELs with those based on 
much fai nter ORLs, has be e n carr ied out for three samples of 
GDPNe jTsamis et al.l [20031 12004 12 PNe, TBLS04 hereafter ; 
ILiu et alj|2004j|bl 12 PNe. LLBL04 hereafter; IWesson et al.ll2005l 
23 PNe, WLB05 hereafter). In addition, case studies for a num- 
ber of individual PNe have been published, some of them show- 
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ing particularly prominent and rich O RL spectra from heavy el 



emen t ions. They i nclude NGC 700 9 dLiu et al.l [l995l: iLuo et al. 1 
200 lh . NGC 6 153 dLiu et alj l2000h . NGC 6543 dWesson & Liul 
2004) , Me 1 -1 dShen et alj2004l) . NGC 150 1 dErcolano et alj2004h 
NGC 7027 dZhang et alj|2005l) and Hf2-2 dLiu et alj|2006h . In the 
current work, six PNe observed belong to the disk population (c.f. 
TableQJ. In all the studies, the nebulae were observed using low- to 
intermediate-resolution long-slit CCD spectrophotometry and then 
analyzed in a similar manner (plasma diagnostics, atomic data and 
ICF scheme). All those disk PNe, including six from the current 
work, are grouped to form a homogeneous sample of GDPNe to 
compare with the current sample of GBPNe. A few nebulae ob- 
served and analyzed using similar techniques were excluded, eithe r 
becau se of their peculiar ities, e.g. Abell30 I Wesson et alj|2003l) , 
Mz 3 dZhang & Liul2002h . M 2-24 dZhang & Liul2003r) . or because 
the objects actually belong to a different population, e.g. NGC 4361 
jLiulll998h and DdDm 1 (WLB05); both are believed to be halo 
PNe. We note that fluxes for a number of lines published by WLB05 
for PN M 3-27 seem to be incorrect. We have therefore remeasured 
the lines and repeated the analysis. In total, our final sample of 
GDPNe consists of 58 PNe and will be referred to hereafter as TLW 
(V for TBLS04, 'L' for LLBL04 and 'W' for WLB05 - the three 
samples constituting the majority of the TLW sample). 

We note that most of the objects observed by TBLS04 are 
within the solar circle, whereas the northern PNe in the WLB05 
sample are mostly outside the solar circle, and the LLBL04 sample 
is a mixture. The mean O/H abundances deduced from the three 
samples are 8.70, 8.58 and 8.60, respectively. The corresponding 
values after excluding Type-I PNe (and the very low excitation PN 
in the WLB05 sample) are 8.74, 8.57 and 8.59, respectively. The 
small variations of mean oxygen abundances between the samples 
thus reflect an abundance gradient across the Galactic disk. We also 
note that the average metal abundances yielded by the KB94 sam- 
ple of disk PNe are systematically higher than those from the TLW 
sample, varying from about 0. 1 dex for oxygen, neon and sulfur to 
about 0.2 dex for carbon, nitrogen and argon. Part of the differences 
may be caused by abundance gradients across the disk. Most PNe 
in the KB94 sample are southern hemisphere nebulae and therefore 
within the solar circle. For the 38 and 16 PNe within and outside 
the solar circle in the KB94 sample, the mean oxygen abundances 
are 8.71 and 8.60, respectively. 



7.1 Comparison between bulge and disk samples 

Table [14] shows that compared to the disk PN sample TLW, the 
average elemental abundances of the bulge PN sample WL07, are 
consistently higher, except for C/H derived from ORLs and Cl/H 
deduced from CELs. If we define AAf(X) as the difference be- 
tween the average abundances of element X relative to H derived 
from disk bulge sample WL07 and from the disk sample TLW, 
AAf (X) = M(X, GBPNe) - Af(X, GDPNe), then AM(X) 
ranges from 0.04 dex for He/H from ORLs, 0. 10 dex for O/H from 
both CELs and ORLs, ~ 0.15 dex for Ne/H and Ar/H from CELs, 
0.15 dex for Mg/H from ORLs, 0.17-0.20 for N/H from CELs and 
from ORLs, and 0.21 dex for S/H from CELs. 

Fig. [21] compares distributions, normalized to unit area, of 
elemental abundances relative to hydrogen and of heavy element 
abundance ratios for the WL07 bulge sample (heavily shaded his- 
tograms) and the TLW disk sample (lightly shaded histograms). 
Panels (a) - (c) of Fig. [21] show that abundances of He and O of 
bulge PNe peak at higher values compared to disk PNe. Panels 
(d) and (g) show that, as already pointed out in £16.21 fc.f . Fig.l 1 9b, 



bulge PNe have lower C/O ratios than disk PNe, by approximately 
0.2 dex. On the other hand, bulge and disk PNe seem to have simi- 
lar distributions of N/O and Ne/O abundance ratios (as discussed in 
%.2\ nitrogen and neon abundances deduced from ORLs and con- 
sequently N/O and Ne/O abundance ratios based on ORLs are less 
reliable compared to those of carbon and oxygen). 

It is generally accepted that abundances of oxygen, neon, sul- 
fur and argon are hardly affected by nucleosynthesis and dredge-up 
processes that occur in progenitor stars of PNe, thus their abun- 
dances should reflect values of the ISM at the time those LIMS 
were formed. For these four elements, their average abundances 
deduced from CELs for the bulge sample WL07 are consistently 
higher by ~ 0.1-0.2 dex than for the disk sample TLW. Corrobo- 
rating evidence is provided by the average O/H abundances derived 
from ORLs for which good atomic data exist and line measure- 
ments are reasonably reliable. The results show that the bulge is 
slightly overabundant with respect to the disk, by approximately 
0. 1 dex. Although the inferred abundance differences between the 
bulge and disk samples are small, the consistency of results for all 
a elements measured as well as from different tracers (CELs versus 
ORLs) strongly suggests that differences are probably meaningful. 

Supporting evidence that the Galactic bulge is more metal-rich 
than the disk has also been found from other independent abun- 
dance studies of bulge PNe. Spectrophotometric abundance anal- 
yses of 30 GBPNe by CMKAS00 show that while the average 
abundances of O, S and Ar of GBPNe are comparable to those of 
GDPNe, there is some evidence that higher elemental abundances 
are slightly more frequently found amongst the bulge PNe. More 
recent work on 44 GBPNe bv lGornv et all d2004l) also came to the 
same conclusion. They compared the average oxygen abundance 
of their sample b (bulge PNe) with sample d (disk PNe within the 
solar circle compiled from the literature) and found that the former 
is about 0.13 dex (or 0.2 dex) more metal abundant than the lat- 
ter. Their findings are roughly in agreement with our results. Note 
that there are only a handful of common objects between the three 
samples of bulge PNe and independent analyses all indicate that 
GBPNe are slightly more metal-rich than GDPNe. 

The lower C/O ratios observed in GBPNe compared to 
GDPNe may suggest that GBPNe evolve from more massive pro- 
genitor stars. As a consequence, a larger fraction of carbon has been 
converted to nitrogen, although it is difficult to rule out selection ef- 
fects as PNe of less massive cores are also less bright and therefore 
harder to observe, particularly those in the distant, heavily extincted 
region of the Galactic bulge. Alternatively, it may indicate that the 
3rd dredge-up process is less efficient for metal-rich GBPNe than 
for GDPNe. In § 17.21 we will discuss in detail the lower C/O ra- 
tios measured for GBPNe and compare the results with some of the 
recent theoretical predictions. 

The distribution of He/H, O/H and Ne/O abundance ratios 
plotted in Fig.|2T]can be well fitted using a single Gaussian profile, 
while those of C/O and N/O appear to be more irregular and double- 
peaked. While observational errors may play a role, the results may 
indicate that there are two populations in each of the two samples. 
A scrutiny of N/H abundances obtained by EBW04 for a sample 
of GBPNe (see the solid histograms in their Fig. 4) also reveals a 
double-peak distribution, peaking at 8.1 and 8.7 dex, respectively, 
similar to what we find for the current sample of GBPNe. Simi- 
larly, N/H abundances of the ECM04 sample also exhibit two local 
maxima, at 7.7 and 8.4 dex, respectively. It is interesting to note 
for each of the three samples of GBPNe, the two local maxima ob- 
served in the N/H abundance distribution are consistent with the 
average N/H abundance ratios for non-Type I and Type-I PNe in 
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Figure 21. Distributions, normalized to unit area, of elemental abundances relative to hydrogen, and of heavy element abundance ratios for GDPNe (lightly 
shaded histograms) and GBPNe (heavily shaded histograms). The ordinate gives the fractional number of nebulae in each bin. The bin size is 0.05 for panel 
(a) and 0. 10 for all the other panels. In each panel, numbers of bulge and disk PNe for which measurements are available are marked from top to bottom. 



the corresponding sample, respectively. The fraction of Type-I PNe 
in the WL07 bulge PN sample is 19 per cent, very similar to the 
fraction of 21 per cent for the TLW disk sample. It seems unlikely 
that the slightly higher metallicity observed for GBPNe compared 
to GDPNe is caused by a mixture of different populations in the 
two samples. 

In a recent theoretical study, Na kasato & Nomotol j2003b have 
modelled the Galaxy formation using three-dimensional, hydrody- 
namic N-body simulations, paying particular attention to the for- 
mation of the bulge component. Their simulations show that the 
bulge consists of two chemically different components - one re- 
sulted from a fast merging of subgalactic clumps in the proto- 
Galaxy while the other formed gradually in the inner disk. Their 



results predict a double-peak [Fe/H] abundance distribution - the 
iron content in the non-merger component is higher than the merger 
component as a result of pollution by Type-la supernovae. It is un- 
clear however how the double-peak [Fe/H] distribution predicted 
for the bulge is related to the double-peak distribution of N/H abun- 
dances observed in GBPNe. 



To summarize, a comparison of elemental abundances of disk 
and bulge PNe shows that GBPNe appear to be more metal-rich 
than GDPNe, by approximately 0. 1 - 0.2 dex. In addition, the mean 
C/O ratio of GBPNe are lower than that of GDPNe, by 0.16 and 
0.27 dex as deduced from CELs and ORLs, respectively. 
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7.2 Carbon and nitrogen 

At the end of the main-sequence evolution when the central hydro- 
gen is exhausted, the star evolves to the red giant branch (RGB). A 
convection zone develops and dredge-ups partially CN-cycle pro- 
cessed material to the surface. This first dredge-up process en- 
hances the surface nitrogen abundance and reduces the surface 
abundance of carbon as well as the 12 C/ 13 C isotope ratio. Fol- 
lowing the exhaustion of central helium and the formation of an 
electron degenerate CO core, the second dredge-up starts in early 
Asymptotic Giant Branch (AGB) stars but only in those more mas- 
sive ones with initial masses in the (uncertain) range from 2.3 to 
8 Mq . The process increases the surface abundances of both he- 
lium and nitrogen and reduces slightly those of carbon and oxy- 
gen. The third dredge-up takes place in AGB stars during the ther- 
mal pulses (helium shell flashes). Helium-burning products, mostly 
12 C and .s-process elements, are brought to the surface. On the way 
to the surface, the material is subject to limited CN-cycle burn- 
ing, allowing for some conversion from C to N and decreasing the 
12 C/ 13 C ratio. In all three episodes of the dredge-up process, the 
surface abundances of He and N are increased, while that of 12 C 
abundance and the 12 C/ 13 C isotope ratio first decrease in the first 
(and second, when it occurs) dredge-up episode and then increase 
in the third. The amount of abundance changes depends mainly on 
the initial mass but also on metallicity of the progenitor star. The 
abundances of C and N relative to other light elements thus provide 
tools to probe the nucleosynthesis and mixing processes as well as 
properties of PN progenitor stars, such as their initial masses and 
metallicities. 

Since the pioneering work bv lRenzini & VoliNl98ll) . the roles 
of LIMS in enhancing the chemical composition of the ISM have 
been studied extensively, both observationally and theoretically. In 
recent work, Marigo (2001, MaOl hereafter) presented quantitative 
predictions of the surface abundance patterns of He, C, N and O for 
different values of initial stellar mass and metallicity and mixing- 
length parameter a. 

In Figs.[22][23]and[24] we plot N/O ratios measured for the 
WL07 sample of GBPNe and for the TLW sample of GDPNe, de- 
duced from CELs and from ORLs against He/H derived from ORL 
analysis and against N/H and O/H obtained from CEL analysis. The 
data are compared with theoretical tracks calculated by MaOl as a 
function of initial mass for initial metallicities Z = 0.04, 0.008 and 
0.019, assuming a = 1.68. 

Fig. [22] shows that the observed distribution of N/O versus 
He/H does not fit well with what is expected, with many data points 
falling to the upper left of the theoretical tracks. Either the theoret- 
ical calculation underestimates the N/O ratio or overestimates the 
He/H enhancement, assuming that most PNe are descendants of 
LIMS, an assumption which is most likely true. The effect is most 
obvious for GDPNe (open symbols). Similarly, Fig. [23] shows that 
while most GBPNe (filled symbols) fall within the band delineated 
by the tracks of various initial metallicities, more than half the open 
symbols fall above and to the left of the track of even the lowest 
metallicity Z = 0.004. The theoretical loci of N/O versus O/H be- 
come almost orthogonal to the abscissa (Fig.l24b for initial masses 
higher than ~ 2 Mq . Again in this diagram, we find that bulge PNe 
are mostly found between the two loci for initial Z = 0.008 and 
0.016, while a significant fraction of disk PNe fall to the left of the 
loci for Z = 0.004. Figs .[22] [23] and [24] thus all seem to indicate 
that bulge PNe originated from environments of higher metallici- 
ties (Z ~ 0.013 from a rough estimate) than disk PNe. In addition, 
Fig.[23] shows that to reach a given value of the N/O ratio, a much 
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Figure 22. N/O abundance ratios derived from CELs (circles) and from 
ORLs (triangles) for the WL07 sample of GBPNe (filled symbols) and the 
TLW sample of GDPNe (open symbols), plotted against helium abundance 
He/H. Overplotted are loci of initial mass predicted by MaOl for initial 
metallicities Z = 0.004, 0.008, 0.019 and mixing length parameter a = 
1.68. The tick marks on the tracks indicate initial stellar mass, ranging from 
~ 0.8 to 5 Mq. Only a few values are labelled. 
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Figure 23. N/O abundance ratios plotted against N/H abundance. For an 
explanation of the different symbols and lines, c.f. caption to Fig. 1221 



higher initial mass is required for a star of higher initial metallicity 
than for one of lower initial metallicity. Fig.1231 shows that almost 
all disk PNe appear to evolve from stars of initial mass lower than 
~ 4 Mq (assuming they have an initial metallicity close to or lower 
than 0.008), while a significant fraction of bulge PNe have initial 
masses higher than this. We thus conclude that PNe of massive pro- 
genitors are more frequently found in the bulge than in the disk, at 
least for the limited sample of bulge PNe analyzed in the current 
work. 

Loci of the C/O ratio versus He/H and O/H abundances cal- 
culated by MaOl are plotted in Figs. [25] and [26] respectively, and 
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compared to observations. For low-mass progenitor stars, in the 
absence of Hot-Bottom-Burning (HBB) which destroys carbon, the 
final surface abundance of carbon as well as the C/O ratio increase 
with increasing initial mass, as a result of the increasing number 
of dredge-up episodes experienced during the Thermal-Pulse-AGB 
phase. The quantities reach maximum values around 2-3 Mq. Be- 
yond that, they start to decline as initial mass increases as a con- 
sequence of less dredge-up events. For even higher initial masses, 
HBB kicks in and eventually prevails. 

As noted earlier ( %.2\ Fig.l 19t. GBPNe appear to have lower 
C/O ratios than GDPNe. Fig.[25]shows that the majority of GBPNe 
have C/O ratios lower than the predicted values of MaOl for initial 
values ranging from ~ 0.8-5 Mq. Even for GDPNe, over half of 
them have C/O ratios lower than the lowest predicted values. Simi- 
lar trends are also found in Fig. [26] Thus, it sees that the theoretical 
calculations of MaOl have overestimated the yield of carbon, by as 
much as about 0.4 dex. The predicted C/O ratios can be lower for 
initial masses higher than 5 Mq . However, it seems highly unlikely 
that over half the PNe in our sample descend from such massive 
stars. 

Fig. [26] can be used to obtain better estimates of the average 
initial metallicity for our disk and bulge PN samples than possible 
with Fig. [24] because the C/O ratio is more sensitive to the initial 
metallicity than N/O is. Fig. [26] shows that essentially all PNe in 
our sample of GBPNe fall between the two loci of Z — 0.008 and 
0.019, while most PNe in the current sample of GDPNe fall close 
to the locus of Z = 0.008. 

In summary, both from the actually observed distributions of 
oxygen abundances (Fig,l21t for bulge and disk PNe as well as from 
the measured N/O and C/O abundance ratios (Figs. I23ll24l andl26ll. 
there is strong evidence indicating that GBPNe were formed in en- 
vironments of higher metallicities (Z ~ 0.013) compared to disk 
PNe. The latter seem to have initial metallicities < 0.008. We also 
find that the current sample of GBPNe may originate from more 
massive stars than the disk PNe. Finally, we point out that MaOl 
may have overestimated carbon yield by over a factor of two. 
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Figure 26. C/O ratios derived from CELs and from ORLs plotted against 
oxygen abundances derived from CELs for GBPNe and GDPNe. See also 
the caption to Fig. l22l for more description. 



7.3 Oxygen, neon and argon 

There is a tight linear relation between neon and oxygen abun- 
dances as illustrated in Fig. [27] The only exception is Sp 4-1 
(WLB05). For this peculiar object, with the exception of oxygen, 
its abundances of all other heavy elements are lower than the aver- 
age of disk PNe by 0.5 to 0.8 dex. Similar correlations between the 
neon and oxygen abundances h ave previously been observed (e.g. 
iHenrvll 19891 iHenrv et al.ll2004T) and interpreted as a direct conse- 
quence of the fact that both neon and oxygen originate from pri- 
mary nucleosynthesis in massive stars (> 10 Mq), and is therefore 
nearly independent of the evolution of LIMS, the progenitor stars 
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Figure 27. Neon abundances plotted against those of oxygen for GBPNe 
(filled circles) and GDPNe (open circles). For both elements, abundances 
deduced from CELs are adopted. The straight line represents a linear re- 
gression for both bulge and disk PNe. 



of PNe. A linear regression of the data points for disk and bulge 
PNe yields Ne/H = (-2.42 ± 0.70) + (1.21 ± 0.08)O/H, with 
a correlation coefficient of 0.85. The result is compatible to what 
Henry (1989) obtained for a sample of Galactic, Magellanic and 
M31 PNe: Ne/H = (-2.14 ± 0.33) + (1 16 ± .04)O/H. Sim- 
ilar results are also obtained by Henry et a for a sample of 
85 Galactic PNe and by ! Stanghellini et alj j200"r3) for a sample of 
79 northern hemisp heric PNe. W e note that from both our linear re- 
gression and that o f lHenryH 19891) for PNe spanning a wide range of 
metallicities, there is marginal evidence that the slope of linear rela- 
tion is larger than unity. To further explore this important issue, we 
have calculated the average Ne/O abundances for SMC and LMC 
PNe, using the data published by llxisv & Denne feld (2006J). For 
the 30 SMC and 120 LMC PNe in their sample, we obtain an av- 
erage Ne/O ratio of 0.154 (or -0.81 dex) and 0.183 (-0.73 dex), 
respectively, significantly lower than the average values of 0.24 
(—0.61 dex) and 0.25 (—0.57 dex) for the Galactic disk and bulge 
samples. The average oxygen abundances for the SMC and LMC 
samples are 8.09 and 8.38, respectively. The results thus strongly 
suggest that there is a time delay of neon enrichment relative to 
oxygen and neon abundance is lower at low metallicities - a result 
that may bear some important implications for our understanding 
of the enrichment of the ISM by massive stars and/or the possible 
evolution of the initial mass function (IMF) with metallicity. 

Similar but looser linear correlations have also been found be- 
tween abundances of oxygen and those of the other two a-elements, 
sulfur and argon. Linear regression yields S/H = (— 2.64±1.32) + 
(1.10 ± 0.15)O/H and Ar/H = (-2.34 ± 1.21) + (0.99 ± 



0.14 )O/H, with a correlation coefficient of about 0.6 in both cases. 

Stanghellin i et al.1 feOOrJ) found an average Ne/O ratio of 
0.27 ± 0.10 for their sample of Galactic PNe. They argued that the 
ratio is independent of the nebular morphology and thus on stellar 
population. For our samples of bulge and disk PNe, we obtain an 
average Ne/O abundance ratio of 0.25 ± 0.08 an d 0.24 ± 0.08, re- 
spectively, in good agreement with the result of [stanghelli ni et all 



(2006). It is interesting t hat the Ne/O abund ance ratios obtained 
for the Orion Nebula by Simp son et alj J 19981) and for two H II 
regions, G333.6-0.2 and W43. by ISimpson et al.1 J2004l) are also 
0.25. The Ne/O abundance ratios derived from PN and H II region 
observat ions are therefore consistently high er than the solar value 
of 0.15 Cpdders 20031: lAsplund et alj|2005l) . If we accept the av- 
erage Ne/O ratio of 0.25 deduced for Galactic PNe, and co mbine 
the va lue with the solar oxygen abundance of 8.69 from Lodders 
we obtain a neon abundance of 8.09, which is significantly 
higher than the currently recommended value of 7.87 for the Sun 
but is remarkably consistent with the value 8.11 deduc ed from solar 
coronal line measurements (Feldman & Widing 2003). 

The relatively high neon ab undance of PNe has also be en ob- 
served and discussed previously. Potta sch & Beintem3 i 19991) sug- 
gested that neon is created in NGC 6302, a bipolar PN ionized 
by a hidden ionizing source of probably the highest temperature 
known for a star (> 200, 000 K). Using multi-waveband data, 
they deduced ionic abundances from Ne + up to Ne 5+ and ob- 
tained a total Ne elemental abundance of 8.34, which is 0.5 and 
0.3 dex higher than the solar valu e and the average abundance of B- 
stars, respectively. Theoretically, Karakas & Lattanzio ( 2003) sug- 
gest that some PNe may have experienced moderate neon enrich- 
ment. The mechanism occurs for stars of initial metallicities less 
than 0.008 and initial masses in the range of 2-4 Mq, where 22 Ne 
is produced in sufficient quantities to affect the total neon abun- 
dance by more than 20 pet cent, about 0.08 dex. However, if the 
mechanism only operates in metal-poor stars in a narrow range of 
initial mass, then it is unlikely it will aff ect the average neon abun- 
dance of PNe by any significant amount. IPottasch & Bernard- Salas 
(2006) argued that the higher Ne/O ratio in PNe compared to the 
Sun is more likely caused by the destru ction of oxygen in progen- 
itor stars of PNe. Model calculations by Karakas (2003) show that 
some oxygen is indeed destroyed in stars more massive than 5 Mq 
of initial metallicity Z — 0.004 and 0.008. The process however 
also produce copious nitrogen, which is inconsistent with observa- 
tions. On the other hand, her calculations for stars of initial metal- 
licity Z — 0.02 (i.e. solar composition) and of initial masses of 6 
and 6.5 Mq, a moderate N/O ratio of about 1.2 is found, which is 
more compatible with observations. To explain the above the aver- 
age neon abundance of PNe requires a significant fraction of PNe 
evolve from progenitor stars in such a restrictive range of initial 
mass and metallicity, an assumption that is obviously unrealistic. 
The average oxygen abundance deduced from CELs for the TLW 
disk sample of PNe is 8.60, whereas that of the WL07 bulge sample 
is 8.70, close to the solar value of 8.69 (c.f. Tablell4t. Thus deple- 
tion of oxygen in PNe is likely to be small if any, and it is unlikely 
that the high Ne/O ratios observed in PNe is caused by the destruc- 
tion of oxygen in their progenitor stars. From the above discussion, 
we conclude that abundances deduced from PNe provide strong, al- 
be it indirect, evidence th at the solar neon abundance recommended 
bv lAsplund et al.1 J2005h may have been underestimated. 

Helioseismological studies have indeed pointed to a higher so- 
lar neon abundance. As a consequence of the application of a time- 
dependent 3D hydrodynamical model of the solar atmosphere in- 
stead of the earlier ID hydrodynamical model, the metal content in 
the solar convection zone has decreased by almost a factor of two. 
Amongst them, abundances of key elements such as, C, N and O 
have been lowered by appr oximately 0.2-0.3 d ex from the earlier 
widely adopted values (see Aspl und et al.l2005l for a review). Those 
large downward revisions however have posed serious challenges, 
most notably for helioseismology. Solar models which adopt the 
new solar abundances as input can no longer satisfactorily match 
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helioseismological observations. On the other hand. iBahcall et all 
show that raising the solar neon abundance to 8.29 (or 
Ne/O = 0.42) can in principle bring the solar model back into 
agreement with the helioseismological measurements. The sugges- 
tion is supported by Chandra measurements of Ne/O ratios for a 
sample of 21 nearby so lar-like stars, yielding an error- weighted 
mean Ne/O ratio of 0.41 jDrake & Test a 2005). However, a reanal- 
ysis of the solar corona data bv lSchmelz et al.ld2005h recover ed the 
old lower solar N e/O ratio of 0.15 (jGre vesse & Sauval 1998j). In a 
more recent studv. lLiefke & Schmitq i 20061) derived a Ne/O ratio of 
~ 0.28, an intermediate value for the inactive, solar-like star a Cen 
(primarily a Cen B, which is the dominant component in X-rays). 
It seems that the problem of the solar Ne/O abundance ratio and 
consequently the absolute abundance of Ne remains to be sorted 
out. 



7.4 Magnesium 

Tablell4lshows that GBPNe have an average CEL abundance of the 
a-element oxygen that is almost identical to the solar value. While 
for the third-row a-element magnesium, its mean abundance de- 
duced from ORLs, 7.71±0.06, is about 0.15 dex higher than the 
solar value 7.55. The latter is almost identical to the average mag- 
nesium abundance 7.56, also derived fro m ORLs, for the TLW disk 
sample of PNe. As previously argued bv lBarlow et al.l j2003l) . de- 
pletion of magnesium onto dust grains in PNe is likely to be small, 
probably less than 30 per cent. 

The result confirms previous studies which find that that 
magnesium is enhanced in the Galactic bulge dMc William & Rich] 
1 1 9941) . Large enhancements of a-element s neon and magnes ium 
are also in the peculiar bulge PN M2-24 fehang & Liull2002l) . In 
a recent abundance study of 27 RGB stars towards the Galactic 
bulge in Baade's Window based on Keck/HIRES echelle observa- 
tions, iFulbright etall d2005l 120061) found enhanced [a/Fe] ratios 
for all bulge giants, including those at super-solar metallicities. The 
enhancement of magnesium in GBPNe and the large [a/Fe] ratios 
in bulge giants, combined with the old age for the bulge, point to 
Type II supernovae as the primary enrichment process in the bulge. 



8 CONCLUSIONS 

In the current work, we present careful abundance analyses for a 
sample of 25 GBPNe plus 6 GDPNe, using both CELs and ORLs. 
The results are compared to those of disk PNe published in the 
recent literature and to the theoretical predictions of MaOl. We find, 

(i) CEL abundances deduced in the current work are consistent, 
within 0.2 dex, with those published in th e literature for samples 
of GB PNe, incl uding W88 dWebsten[l988l) . RPDM 97 JiTatag et al ' 



1997h . SRM98 JStasinska et aill 998). CMKAS00 dCuisinier et al 



20001) , ECM04 l lEscuderoetalJ|2004l) and EBW04 fexteretal 
2004). An exception is EC01 jEscudero & Cost A 1200 lh . which 
yields average abundances that are at odds with all other samples 
of GBPNe. 

(ii) GBPNe are found to be slightly more metal rich on average 
than GDPNe, by approximately 0.1 to 0.2 dex. GBPNe also have 
C/O ratios systematically lower than GDPNe, by approximately 
0.2 dex. 

(iii) Bulge PNe, similar to disk PNe, also exhibit ubiquitous 
discrepancies between T c and heavy element abundances derived 
from CELs and those deduced from CELs. The distribution of the 



abundance discrepancy factor, adf, for the bulge PNe are similar to 
that of disk PNe. Analyses also confirm many of the trends previ- 
ously found for disk PNe, such as the adf is positively correlated 
with the difference between the [O III] forbidden line temperatures 
and the H I Balmer jump temperatures, and that for a given neb- 
ula, values of adf deduced for all abundant second-row elements 
are comparable. 

(iv) Based on ORL abundance analyses, we have obtained reli- 
able C/O abundance ratios for a substantial sample of GBPNe. By 
comparing the observed C/O ratios with the recent theoretical pre- 
dictions of MaOl, we conclude that GBPNe were probably formed 
in more enriched environments than disk PNe, and that MaOl may 
have overestimated the carbon yield by as much as 0.4 dex. In ad- 
dition, we find that GBPNe tend to have more massive progenitor 
stars than GDPNe do. 

(v) The average Ne/O ratio deduced from PNe in the WL07 and 
TLW samples is about 0.25, much higher than the current recom- 
mended solar value. This may points to a solar neon abundance that 
has been underestimated by 0.2 dex. 

(vi) For the first time, abundances of the important a-element 
magnesium have been determined using ORLs for a handful of 
bulge and disk PNe. We find that GBPNe have an average abun- 
dance about 0.13 dex higher than GDPNe. The latter have an aver- 
age abundance almost identical to the solar value. The enhancement 
of magnesium in GBPNe and [a/Fe] in bulge giants suggests that 
the primary enrichment process of the bulge was Type-II SNe. 
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APPENDIX A: LINE FLUXES 



Table Al. Ultraviolet and optical line fluxes on a scale where H/3 = 100. 
A superscript s denotes the emission is mainly from the central star and a 
ps indicates that stellar emission contributes at least partially the observed 
feature. (This is only a sample of this table: the full version is avaiable 
online as supplementary meterial.) 



Cn 1-5 



^obs 


F(A) 


/(A) 


Ion 


Ao 


Mult 


1246.53 s 


11.480 


72.136 


N V 


1240 




1345.60 s 


4.309 


20.958 


O VI 


1342 




1375.11 s 


8.256 


37.868 


O V 


1371.29 




1553.86?" 


56.230 


213.618 


C IV 


1551 




1641. 53 ps 


16.560 


59.797 


Hell 


1640 




1720.70 s 


6.039 


21.368 


NIV 


1718.55 




1908.64 


36.220 


145.257 


C III] 


1906.68 




3189.71 


3.944 


6.356 


He I 


3188.67 


V3 


3205.03 s 


2.797 


4.472 


HeH 


3203.10 


3.5 


3704.33 


1.703 


2.284 


H 16 


3703.86 


H16 


3712.44 


1.338 


1.792 


H 15 


3711.97 


H15 


3722.31 


20.510 


27.409 


H 14 


3721.94 


H14 


3727.06 


58.640 


78.365 


[Oil] 


3726.03 


Fl 


3734.84 


1.747 


2.329 


H 13 


3734.37 


H13 


3750.62 


1.950 


2.594 


H 12 


3750.15 


H12 


3770.04 


2.852 


3.777 


H 11 


3770.63 


Hll 


3799.11 


5.574 


7.341 


H 10 


3797.90 


H10 


3817.55 


4.134 


5.420 


He I 


3819.62 


V22 


3833.94 


8.317 


10.867 


H9 


3835.39 


H9 


3867.80 


70.260 


91.077 


[Ne III] 


3868.75 


Fl 


3887.79 


18.000 


23.228 


H8 


3889.05 


H8 


3966.49 


23.230 


29.474 


[Ne III] 


3967.46 


Fl 


3968.79 


11.940 


15.132 


H7 


3970.07 


H7 


4008.82 


0.236 


0.296 


He I 


4009.26 


V55 


4068.16 


3.894 


4.820 


[sm 


4068.60 


Fl 


4071.56 


0.204 


0.252 


O II 


4072.16 


V10 


4075.82 


1.383 


1.708 


[S II] 


4076.35 


Fl 


4088.94 


0.055 


0.068 


on 


4089.29 


V48a 


4096.90 


0.559 


0.687 


NIII 


4097.33 


VI 


4101.27 


20.530 


25.210 


H6 


4101.74 


H6 


4120.15 


0.292 


0.357 


Oil 


4121.46 


V19 


4132.42 


0.087 


0.106 


Oil 


4132.80 


V19 


4143.32 


0.449 


0.546 


He I 


4143.76 


V53 


4152.80 


0.100 


0.121 


O II 


4153.30 


V19 


4155.96 


0.068 


0.082 


Oil 


4156.53 


V19 


4236.61 


0.049 


0.058 


N II 


4236.91 


V48a 


4241.25 


0.113 


0.133 


N II 


4241.24 


V48a 


4266.66 


1.109 


1.305 


C II 


4267.15 


V6 


4275.85 


0.092 


0.108 


Oil 


4275.55 


V67a 


4283.94 


0.049 


0.058 


O II 


4283.73 


V67c 


4339.93 


40.520 


46.763 


H5 


4340.47 


H5 


4349.05 


0.123 


0.142 


Oil 


4349.43 


V2 


4362.69 


3.061 


3.509 


[O III] 


4363.21 


F2 


4387.39 


0.719 


0.819 


He I 


4387.93 


V51 


4391.59 


0.043 


0.049 


Ne II 


4391.99 


V55e 


4437.13 


0.066 


0.074 


He I 


4437.55 


V50 


4470.96 


5.974 


6.650 


He I 


4471.50 


V14 


4480.86 


0.021 


0.023 


Mg II 


4481.21 


V4 


4530.00 


0.053 


0.058 


NIII 


4530.86 


V3 


4561.97 


0.076 


0.082 


Mg I] 


4562.60 




4570.59 


0.221 


0.240 


Mg I] 


4571.10 




4630.08 


0.134 


0.143 


Nil 


4630.54 


V5 


4633.66 


0.294 


0.313 


NIII 


4634.14 


V2 


4639.80 


0.457 


0.486 


NIII 


4640.64 


V2 


4641.01 


0.540 


0.574 


NIII 


4641.84 


V2 


4648.86 


0.378 


0.401 


Oil 


4649.13 


VI 


4657.71 


0.941 


0.995 


[Fe III] 


4658.10 


F3 


4661.39 


0.081 


0.086 


Oil 


4661.63 


VI 
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Table A2. Infrared line fluxes of PN Cn 1-5, M 2-23 and NGC 6567. 



Lines 


F{\) 


7(A) 




(10" 12 ergscm- 2 s- 1 ) 


[/(H/3) = 100] 


Cn 1-5 






[O III] 52/im 


21.16 


343 


[O III] 88/im 


8.39 


136 


[N III] 57/im 


9.66 


157 


[C II] 158/im 


0.75 


12 


M2-23 






Hi 6-4 2.63^m 


0.75 


285 


Hl5-4 4.05/im 


1.37 


520 


Hi 6-5 7.46/im 


0.45 


171 


[Arlll]8.99/im 


1.06 


401 


[S iv]10.5;um 


2.43 


923 


[Nell]12.8/im 


0.79 


299 


[Nelll] 15.6/im 


13.8 


5241 


[S III] 18.7//m 


0.85 


324 


[Nem]36.0/im 


0.64 


243 


NGC 6567 






Hl5-4 4.1^m 


0.81 


7.07 


[Arlll]8.99/im 


1.08 


9.39 


[S IV] 10.5/um 


10.1 


87.7 


[Nelll] 15.6/j.m 


28.9 


252 


[S III] 18.7/im 


1.25 


10.85 


[Nem]36.0/im 


2.25 


19.6 


[Om] 52//m 


51.1 


455 


[O III] 88/zm 


29.7 


264 



